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Abstract. We study the logarithmic and ratio asymptotic of linear forms constructed from 
a Nikishin system which satisfy orthogonality conditions with respect to a system of measures 
generated from another Nikishin system. This construction combines type I and type II multiple 
orthogonal polynomials. The logarithmic asymptotic of the linear forms is expressed in terms 
of the extremal solution of an associated vector valued equilibrium problem for the logarithmic 
potential. The ratio asymptotic is described by means of a conformal representation of an 
appropriate Riemann surface of genus zero onto the extended complex plane. 
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Let s be a finite positive Borel measure supported on a compact subset supp(s) of the real 
line, and (wq, ... , k;* ), (u)q, tujL,) ^ e t w0 systems of continuous functions on supp(s). Fix 
ni = (m >0 , ni,i, ni )OTl ) e Z™ 1+1 and n 2 = (n 2 ,o, n 2 ,i, ■ ■ ■ , n 2 , m2 ) S Z™ 2+1 . Set |m| = 
n i,o + n i,i + ■ • • + Hi, mi , |ri2 1 = + • • ■ + «2,m 2 ; an d n = (ni;n2). In the sequel, we suppose 
that |n 2 | + 1 = |ni|. 

Let | ni | > 1- It is easy to see that there exist polynomials a n .o, a n ,i, ■ ■ ■ , a„. mi such that: 

i) deg(a nj ) < riij — 1, j = 0, . . . , mi, not all identically equal to zero. 

ii) For k = 0, . . . , m 2 



(deg(a nj -) < —1 means that a nj = 0.) 

When m 2 = the polynomials (a ni o, ■ • ■ ,a n . mi ) are called type I multiple orthogonal polyno- 
mials. If mi = 0, a nj o is called a type II multiple orthogonal polynomial. The case mi = m 2 = 

The first three authors in alphabetical order were supported by grants MTM 2006-13000-C03-02 of Ministerio 
de Ciencia y Tecnologi'a and CCG07-UC3M/ESP-3339 of Comunidad Autonoma de Madrid-Universidad Carlos III 
de Madrid. V. N. Sorokin received support from grants RFBR-08-01-00317 and NSh-3906. 2008.1. 



A. M.S. Subject Classification. Primary: 30E10, 42C05; Secondary: 41A20. 



1. Introduction 





2 U. FIDALGO, A. LOPEZ, G. LOPEZ, AND V. N. SOROKIN 

reduces to the usual definition of orthogonal polynomial. When mi, rri2 > 1 these multiple orthog- 
onal polynomials are called of mixed type. 

Multiple orthogonal polynomials appear in problems connected with the algebraic independence 
of functions and numbers (type I) and in questions related with simultaneous rational approxima- 
tion (type II). Those of type II are formed by polynomials which share orthogonality conditions 
with a system of measures which may be written in the form of orthogonality relations with respect 
to a family of generalized polynomials (that in the sequel we call linear forms). In type I, the linear 
forms arc defined through full orthogonality relations with respect to a single measure. Mixed type 
multiple orthogonal polynomials occur in stochastic models connected with random matrices and 
non intersecting random paths, see [6]. Mixed type multiple orthogonal polynomials as presented 
above were considered in [24] and their algebraic properties studied in |26] . 

We will restrict our attention to mixed type multiple orthogonal polynomials in which the 
linear forms are generated by two (not necessarily distinct) Nikishin systems of measures. Nikishin 
systems of measures were introduced in [18] . Before going into details let us mention some papers 
which constitute our starting point. 

E. M. Nikishin studied the asymptotic behavior of the linear forms generated by a Nikishin 
system of measures in [19] (see also [15] and the last section in [20]). He described the logarithmic 
asymptotic of type I multiple orthogonal polynomials in terms of the solution of a vector equilibrium 
problem for the logarithmic potential. Later, Gonchar-Rakhmanov-Sorokin studied in [TT] the rate 
of convergence of Hermite-Pade approximation of generalized Nikishin systems of functions and the 
logarithmic asymptotic of their associated type II multiple orthogonal polynomials. The solution is 
also characterized by a similar vector equilibrium problem. In |25j . V. N. Sorokin defines mixed type 
multiple orthogonal polynomials for two Nikishin systems and gives their logarithmic asymptotic. 

Let s be a finite positive Borel measure supported on a bounded interval A of the real line K 
such that s' > almost everywhere on A and let {Q n },n € Z+, be the corresponding sequence 
of monic orthogonal polynomials; that is, with leading coefficients equal to one. In a series of two 
papers (see [21] and [22]), E. A. Rakhmanov proved that under these conditions 

(1) ito%*M = J^) £ CC \A 

«ez+ Q n {z) (p'(oo) 

(uniformly on each compact subset of C \ A), where ip(z) denotes the conformal representation 
of C \ A onto {w : \w\ > 1} such that <p(oo) = oo and y'(oo) > 0. This result attracted great 
attention because of its theoretical interest within the general theory of orthogonal polynomials 
and its applications to the theory of rational approximation of analytic functions. Simplified proofs 
of Rakhmanov's theorem may be found in [53] and [TB] • 

This result has been extended in several directions. Orthogonal polynomials with respect to 
varying measures (depending on the degree of the polynomial) arise in the study of multipoint Pade 
approximation of Markov functions. In this context, in [12) and [13j . an analogue of Rakhmanov's 
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theorem for such sequences of orthogonal polynomials was proved. Recently, S. A. Denisov [7] 
(see also [17]) extended Rakhmanov's result to the case when supp(s) = A U e C 1, where A is a 
bounded interval, e is a set without accumulation points in R\ A, and s' > a.e. on A. A version 
for orthogonal polynomials with respect to varying Denisov type measures was given in [2]. 

For multiple orthogonal polynomials associated with Nikishin systems of measures an analogue 
of Rakhmanov's theorem was proved in [1] and extended in [14] to the case when the measures in 
the Nikishin system are as those considered by Denisov. 

Let us define the notion of Nikishin system of measures. Let a a , erg be two measures with con- 
stant sign supported on R and let A Q , A^ denote the smallest intervals containing their supports, 
supp(<7 a ) and supp(erg), respectively. We write Co(supp(cr Q )) = A a . Assume that A a n A^ = 
and define 

(cr Q , (jp){x) := J -^r^j d<J a (x) = ap(x)da a (x). 

Therefore, (a a , ap) is a measure with constant sign and support equal to that of a a . 

For a system of intervals Ao, . . . , A m contained in R satisfying Aj nA ]+ i = 0, j = 0, . . . , to — 1, 
and finite Borcl measures ao, ■ • • , c m with constant sign in Co(supp(<7j)) = Aj, such that each one 
has infinitely many points in its support, we define recursively 

(<7 , (Tl, . . . , (Tj) = (cr 0) (<ri, . . . , aj)), j = 1, . . . , to. 

We say that (so, . . . , s m ) — N(ao, . . . , cr m ), where 

so = (o"o) = °o, si — (ao, ax), ■ ■ ■ , s m = (a , . . . , a m ) 

is the Nikishin system of measures generated by (ao, . . . , a m ). In the sequel, when referring to a 
Nikishin system the condition Aj n A^+i = 0,j = 0, . . . ,m — 1, is always assumed to hold. Notice 
that all the measures in a Nikishin system have the same support, namely supp(cro). We will 
denote (sjj = aj) 

Sj,k = (cj, • • • j < j < k < to. 

Take two systems S 1 = (sJ,...,s^J = Af(a^ , . . . , a^ ), S 2 = (sg, . . . , s^J = N(al , . . . , a^ 2 ) 
generated by mi + 1 and TO2 + 1 measures, respectively. The two systems need not coincide, but we 
will assume that (Tq = a^; that is, both systems stem from the same basis measure. The smallest 
interval containing supp(<7*) will be denoted Co(supp(crp) = A*. 

Fix ni = (m,o, ni t x, n\, mx ) € Z" ll+1 and n 2 = (n 2j0 , n 2 ,i, • ■ • , n 2 , m2 ) € Z" l2+1 . Set |ni| = 
^1,0 + + • • • + ni,mi) |n 2 1 = n 2 ,o + • • • + n 2jm2 , and n = (m; n 2 ). In the sequel, we suppose 
that |n 2 | + 1 = |nxj. 

Let | ni| > 1. The system of polynomials a nj o, a ni i, . . . , a niTni satisfying: 

i') deg(a nj ) < nij — 1, j = 0, .. . , mi, not all identically equal to zero. 
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ii') For k = 0, . . . , mi 

(2) J x u fa n>0 (a:) + a a ^{x)s\^{x) j dsl k (x) = 0, v = 0, . . . , n 2 ,k - 1, 

(deg(a nj j) < —1 means that a n j = 0) is called a system of mixed type multiple orthogonal 
polynomials relative to the multi-index n = (ni;n2) and the pair (5' 1 ,5 2 ) of Nikishin systems. 
This concept was first introduced in [2"5] . 

Finding a n ^, . . . , a n mi reduces to solving a homogeneous linear system of |n 2 | equations on |n! | 
unknowns. Since |n.2 1 = | ni | — 1 a non-trivial solution is guaranteed. 

A multi-index n — (ni;ri2) is said to be normal if every solution to i')-fi') satisfies dega nj = 
n\ t j — 1, j = 0, . . . , m. If n is normal, it is easy to verify that the vector (a n ,o, • ■ ■ , dn,mi) is uniquely 
determined except for a constant factor, and in that case we normalize it to be "monic" meaning 
by this that its last entry different from zero has leading coefficient equal to 1. Set 

z™ 1+1 (.) = {m e z™ 1+1 : nx, > • • • > m, mi }. 

In Proposition [U we prove that all n = (ni;n.2) € Z" 11+1 (») x Z™ 2+1 («) are normal. For the 
sequences of multi-indices we shall consider, for almost all n we will have that n\ >mi > 1 and a 
"monic" (<z n ,0) a n,ij ■ • ■ , a n,mi) w ih have a n mi monic. 

Theorem 1 gives the rate of convergence of the |ni|-th root of the linear forms 

mi 

An,o(z) = a n . (z) + ^ a n,k( z )s\,k( Z ^ 
k=l 

under mild conditions on the sequence of multi-indices and the measures generating both Nikishin 
systems. A measure a is said to be regular if 

lim n\l n = l/cap(supp(er)), 

n — >oo 

where cap(-) denotes the logarithmic capacity of the Borel set (•) and K n denotes the leading 
coefficient of the nth orthonormal polynomial with respect to a. For different equivalent forms of 
defining regular measures see sections 3.1 to 3.3 in [28] (in particular Theorem 3.1.1). For short, 
we write (S 1 , S 2 ) S Reg to mean that all the measures which generate both Nikishin systems are 
regular and their supports are regular compact sets. Recall that a compact set is regular when the 
Green's function with singularity at oo of the unbounded connected component of the complement 
of the compact set can be extended continuously to all C. Before stating Theorem 1, we need to 
introduce some notation and results from potential theory. 

Let Ek, k — —m2, ■ ■ ■ be (not necessarily distinct) compact subsets of the real line and 
C = (cj.fe), — W2 < j,k < mi, a real, positive definite, symmetric matrix of order mi + m% + 1. C 
will be called the interaction matrix. By JA(Ek) we denote the class of all finite, positive, Borel 
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measures with compact support consisting of an infinite set of points contained in and M. \ (Ek ) 
is the subclass of probability measures of M(Ek). Set 

Mi = M!{E- m2 ) x • ■ ■ x M!(E mi ). 

Given a vector measure [i = (/i_ m2 ,..., /z mi ) £ Mi and j = -m2,...,mi, we define the 
combined potential 

mi 

(3) W?(x)= ]T c ^ h (^), 

k — — 7U2 

where 

V^(x) = J logj-l—dfikit), 
denotes the standard logarithmic potential of fik ■ We denote 

u% = inf{ Wf{x) :x£Ej}, j = -m 2 , . . . , m x . 

In Chapter 5 of |20j the authors prove (we state the result in a form convenient for our purpose). 

Lemma 1. Assume that the compact sets Ek,k — — m2,...,mi, are regular with respect to the 
Dirichlet problem. Let C be a real, positive definite, symmetric matrix of order mi + mi + 1. // 
there exists Jl — (jl_ m2 , . . . ,JI mi ) £ Mi such that for each j = —m 2 , • ■ ■ , mi 

Wj"(x) = u>j , x £ supp(7* i j) , 

then Jl is unique. Moreover, if Cj,k > when Ej n Ek ^ 0, then Jl exists. 

For details on how this lemma is derived from [301 Chapter 5] see [3J Section 4]. The vector 
measure Jl £ Mi is called the equilibrium solution for the vector potential problem determined by 
the interaction matrix C on the system of compact sets Ej , j = — mi, . . . , m\ . 

Let A = A(pi. o, . . . , pi,m 1 > P2.0, ■ ■ ■ , Pi,m 2 ) C Z™ 1+1 (») x Z™ 2+1 (») be an infinite sequence of 
distinct multi-indices such that 

^ T~H = V\,i € (0,1), j = 0, ...,mi, lim r^- =p 2 ,j e (0,1), i = 0,...,m 2 . 

neA |m| neA |m| 

Obviously, pi, > • • • > Pi, mi ,P2,o > • • ■ > P2,m 2 , and X^qPij = J2T=oP 2 d = 1 - Set 



mi m2 
= X! Pl ' fe ' J = 0, . . . ,T7T-l, P_j = ^P2,fe, j = 0, . . . ,m 2 . 
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Let us define the interaction matrix C which is relevant for the rest of the paper. Take the 
tri-diagonal matrix 

/ p2 ^ _ P- m2 P- m2+1 _ ^ 



(4) 



pi 

— 771; 



-2 + 1 





P-m 2 + lP-™ 2 + 2 



m 2 + l 



-2+ 2 



P 5 



m 2 +2 



V o o o ••• p^j 

This matrix satisfies all the assumptions of Lemma [T] on the compact sets Ej — supp = 
0, 1, . . . , mi, .Ej = supp((T^ J ), j = 0, — 1, . . . , — 77^2 , including > when £j D -Efc 7^ (recall that 
°o = °o)> an0 - it is positive definite because the principal section C r ,r — 1, . . . , mi + 777,2 + 1, of C 
satisfies 



(1 -k 



det(C r ) = p2 ...P^ _ ld et 




\ 



\ 





> 0. 



/ 



Let fJ-(C) be the equilibrium solution for the corresponding vector potential problem. We have 

Theorem 1. Let A = A(pi ;0 , ■ ■ ■ ,Pi, mi ',P2,o, ■ ■ ■ >P2,m 2 ) C Z™ 1+1 (.) x Z™ 2+1 (.), (5 1 , S 2 ) G Reg, 
S* 1 = jV(<7q, . . . , c^i)' arl ^ ^ 2 = ■^'( <7 0) • • • j <T m 2 ) ^ e gi ven - Then 

(5) lim |^ n . (z)| 1/|ni1 = G (z), /CcC\(AjuA}), 

uniformly on each compact subset IC C C \ (Aj U A}), where 

G (z) = exp ( Pi^(z) - V^(z) - 2J2 ^ ) • 



fc=i 



ju = /7(C) = (J2_ m2 , ■ . • ,~p mi ) is the equilibrium vector measure and {ojt m , . . . , w^, ) is f/ie system 
of equilibrium constants for the vector potential problem determined by the interaction matrix C 
defined in Q on the system of compact sets Ej — supp(crj), j = 0, . . . ,m\,Ej — supp(cr 2 _ J ), j = 
-77l 2 , ... ,0. 



Throughout the paper, the notation 



\\mg n (z) = g(z), 

riGA 



KcO, 



stands for uniform convergence of the sequence {gn},n € A, to g on each compact subset IC 
contained in the indicated region (in this case f2). 
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For the next result, we assume that supp(erj) = A* Ue*-, i = 1, 2, where AJ is a bounded interval 
of the real line, > a.e. on A*-, and e*- is at most a denumerable set without accumulation 

points in R\ Aj. We denote this writing S 1 = jV'(cr^ . .. , a} ni ),S 2 = J\f'(a%, . . . , cr^J. Fix a 
vector £ := (Zi;^) where < l\ < mi and < Z 2 < to 2 . We define the multi-index n' := 
(ni + e (l ;ri2 + e' 2 ) = (n^jnj ), where e'* denotes the unit vector of length m, + 1 with all 
components equal to zero except the component {U + 1) which equals 1. It is always assumed that 
both n and n* belong to Z" ll+1 (.) x Z™ 2+1 (.). 

Theorem 2. Let S 1 = Af'(a^, a l mi ), S 2 = N'{al, a 2 n J, and A C Z™ 1+1 (.) x Z™ 2+1 (.) is 

an infinite sequence of distinct multi-indices such that 

sup((m 2 + l)n-2,o - |n 2 |) < oo, sup((mi + l)ni,o - |ni|) < oo. 

n£A nSA 

Assume that there exists I — (Zi ; Z2), < l\ < mi,0 < I2 < m,2, such that for all n € A we have 
that n' = (nj 1 ; n 2 2 ) G Z™ 1+1 (.) x Z™ 2+1 (.). T/ien 

(6) lim = ^O*), £ C C\ (supp(^) U suppK 1 )) , 

neA A n fi{Z) 

where Aq (z) is a one to one analytic function in C \ (Aj U Aj). 

An expression for J${z) will be given in Theorem [7] The answer depends of the conformal 
representation of an associated Riemann surface with m\ + m 2 + 2 sheets and genus zero onto the 
extended complex plane. The previous result is already new when m 2 = and m\ > 2. 

Besides normality, in Section 2 we obtain the orthogonality relations satisfied by the linear forms 
involved in the construction. Section [4] is devoted to the study of the asymptotic distribution of 
zeros of a system of linear forms associated with -A n ,o that allows to prove Theorem[5]in Section[5]of 
which Theorem [T] is a corollary. Theorem 5 was first stated in [24, under the stronger assumptions 
S 1 = N'(o-q, . . . , cr^J, S 2 = M'(o- 2 , . . . , o~ 2 l2 ), and supp(crj) = A*, i = 1, 2; and the proof was only 
carried out for m\ = m 2 = 1. 

In Section [3] we study the interlacing properties of the zeros of the linear forms which is needed 
for the proof of Theorem [7] in Section [6] from which Theorem [2] follows. Section 7 contains a Markov 
type theorem for mixed type Hermite-Pade approximation and some reinterpretation of the theory 
developed in the context of systems of bi-orthogonal linear forms. 

2. Normality and orthogonality relations 

Recall that 

8j,k = {o~j, ■ ■ • , 0-fc), < j < k < m, s 3ij = o-j . 
We denote (s} +l j (z) = l,^ n>mi = a„, mi ) 

n 1 1 

A n ,j(z) := ^a n , fc (z)s) +lfc (z), j = 0,..., mi. 
k=j 
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In [18], E. M. Nikishin introduced the following definition. 

Definition 1. A set of real continuous functions uq(x), . . . , u mi {x) defined on an interval A, is 
called an AT-system for TH\ — (nifi, . . . , n\, mi ) € Z" ll+1 , if for any polynomials ho ■ ■ ■ , h mi such 
that deg(hi) < n\ i — 1, i = 0, . . . , m\, not simultaneously identically equal to zero, the function 

h {x)u (x) H h h mi (x)u mi (x), 

has at most \rii\ — 1 zeros on A (deg(hj) < —1 means that hj = 0). 

Let Z™ 1+1 (*) be the set of multi-indices given by 

Z™ 1+1 (*) = { ni e Z™ 1+1 : fii<k <j such that n, < nj < n k }. 
In connection with AT-systems, in [8] U. Fidalgo and G. Lopez proved 

Lemma 2. Let ni G Z" ll+1 (*) and (s\, . . . , s rni ) = M{<J\, . . . , cr mi ), then the system of functions 
(1, si, . . . , s mi ) defines an AT-system with respect to n\ = (n^o, ■ ■ ■ , Hi jm J on an?/ interval disjoint 
from Co(supp(o"i)). 

For each j = 0,...,mi - 1, we have that (s} +lj+1 , . . . ,s] +l mi ) = N{a) +1 , . . Using 
Lemma [2] it follows that for ni £ Z™ 1+1 (») C Z" ll+1 (*) the linear form A n j cannot have more 
that Nx t j — 1 zeros on any interval disjoint from Aj +1 , where 

Nxj = N ld (n) = mj H h n.i, mi . 

Obviously, the same is true for the polynomial -An. mi = &n,mi • 

Below, we also use the previous 

lemma for linear forms generated by the second Nikishin system. 

Notice that dsy k (x) = ~s\ k (x)do-Q(x). On the other hand, we can replace x v by any polynomial 
of degree < n-2,fc — 1 inside the integral in ([2]). Set 

m 2 

g n 2 (z) = }] K 2 ,k( z )si.k( z )> deg6 n2i /c < n 2 , fc - 1, fc = 0, ...,m 2 . 

fc=0 

(sf (z) = 1). Then © is equivalent to 

(7) J B B2 {x)A n ,o(x)dal(x) = 0, 

for all B n2 as indicated. 

Suppose that *4 n ,o has at most N < |ni| — 1 = |ri2 1 sign changes on the interval Aq. Choose 
polynomials 6 n2 ,fc conveniently so that B n2 changes sign exactly at those points where *4 n ,o changes 
sign on Aq and has a zero of order |ri2 1 — 1 — N at one of the extreme points of A J = Aq. By 
Lemma [2 the linear form B n2 has on Aq at most |ri2 1 — 1 zeros, thus it can only have those zeros 
which we have assigned to it. The continuous function B n2 A n fl has constant sign on Aq. This 
contradicts (JT]). 
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We have proved that A n ,o has |ni| — 1 zeros with odd multiplicity in the interior of Aq = Aj. 
In connection with intervals of the real line, the interior refers to the Euclidean topology of K. In 
short, we shall see that -4 n ,o has no other zeros in C \ A} and that they are all simple. Before 
proving this, let us turn to the question of normality. 

Proposition 1. Let n £ Z™ 1+1 (.) x Z™ 2+1 (.), S 1 =N'(a^,..., a l mi ), and S 2 = Af(a 2 , . . . , a^J, 
be given. Then, n is normal and (a nj o, ■ ■ • , a„ im J is uniquely determined except for a constant 
factor. 

Proof. Assume that there exists j £ {0, . . . , mi} such that dega nj < n\j — 2. Then ni — e J £ 
Z™ 1+1 (*), where e J ' denotes the mi + 1 dimensional unit vector with all components equal to zero 
except the component j + 1 which equals 1. According to Lemma [2] applied to ni — e J , the linear 
form A n .o has at most |ni | — 2 zeros on Aq, but we pointed out before that it has at least |ni | — 1 sign 
changes on this interval. This contradiction yields that for all j £ {0, . . . , mi}, dega nj = — 1, 
which implies normality. 

Now, let us assume that (a nj o, ■ ■ ■ , a n ,mi) and (a* , ■ ■ ■ , a* m ) solve i' and these vectors 
are not collinear. According to what we just proved, for all j £ {0, . . . , mi}, dega nj = dega* j = 
m,j — 1. Take A £ C \ {0} such that deg(a n .o — Aa* ) < Ui,o ~ 2. Obviously, the vector (a n .o — 
Aa* q , . . . , a n . mi — Aa* ) is not identically equal to zero and also solves i ' )— ii' ) which is not possible 
since all non trivial solutions must have all components of maximal degree. □ 

Proposition[T] allows us to determine the "monic" (a nj o, a n ,i, ■ ■ ■ , ffln.mj uniquely and we impose 
this normalization . The next lemma will be used on several occasions. 

Lemma 3. Let Sk,k — l,...,m, be finite signed Borel measures, Co(supp(sfc)) =Acl. Let 
F{z) = fo{z) + J2k=i fk( z )sk(z) £ H(C \ A), where f k £ H(V), k = 0, . . . ,m, and V is a neigh- 
borhood of A. IfF(z) = 0{l/z 2 ),z -> oo, then 

(8) / fk(x)ds k (x) = 

k=i^ 

and F{z) = 0(l/z), z — > oo, implies that 

fk{x)ds k {x) 



III- n 

(9) F(z) = Y,j 



z — x 



k=l ■ 

■2\ 



Proof. Let r C V be a closed smooth Jordan curve that surrounds A. If F(z) — 0(1/ z ), z — > oo, 
from Cauchy's theorem, Fubini's theorem and Cauchy's integral formula, it follows that 

n Til n TfL n n n / \ 7 TTb n 

0= / F(z)dz = V / f k (z)s k (z)dz = y i / / ^^ds k (x) = 2tu V / f k (x)ds k (x), 
Jr k=i Jr k=i j Jt z ~ x k=i J 

and we obtain ([5]). On the other hand, if F(z) — 0(1/ z), z — > oo, and we assume that z is in the 

unbouded connected component of the complement of T, Cauchy's integral formula and Fubini's 
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theorem render 



= J_ f F(()d( = J_^f fk(C)M()d( 



2wi J T z — ( 2iri ^— J J T z — ( 

fk(()d( . vW fk(x)ds k (x) 



k=l" 1 v y fc=l' 

which is ©. □ 

In the sequel, n = (m; n 2 ) € Z" ll+1 (») x Z™ 2+1 (») and m, mi > 1. For j = 0,... 1 m 1 , let Q nj 
be the monic polynomial whose zeros are those of the linear form A n j in the region C \ Aj +1 , 
counting multiplicities (A^ i+1 = 0). In particular, A n ^ ni = a nmi = Q n . mi . From the previous 
proposition, if ?ii, mi > 1, oo is not a zero of any one of these linear forms; thus, oo cannot be 
an accumulation point of such zeros. Though it is not the case, in principle, some of these linear 
forms may have an infinite number of zeros which accumulate on the boundary of the corresponding 
region of meromorphy. In that case, for the time being, Q n ,j denotes a formal infinite product. 

The next proposition is adapted from [15] . 

Proposition 2. Let n e Z™ 1+1 (.) x Z™ 3+1 («), m, mi > 1, S 1 = jV(c^,..., ^J, and S 2 = 
M{oq, . . . , o~ 2 n2 ) be given [recall that Cq = a 2 ). Then, degQ nj = N\,j — l,i = 0, 1 . . ., mi, iis 
zeros are simple and lie in the interior of Aj . If I denotes the closure of a connected component 
of \ supp(erj) then Q n j has at most one zero in I. Moreover, 
r da) (x) 

(10) / x v A aJ (x) 3 \' =0, v = 0,...,N ld -2, j = l,..., mi , 

and for any polynomial q, degq < Nij+i — 1, 

g(A^M = r g(x)A a , j+1 (x) fo^x) j = mi _ h 

Qn,j(z) J Qnj(x) Z-X 

Proof. Using induction on j, we will prove simultaneously the general statement concerning the 
zeros and (|10|) . Then, we prove that on any interval I there is at most one zero of Q n .j- Finally, 
we obtain (fTTj) . For j = 0, we already proved that A n ,o has Ni t o — 1 = |ni| — 1 sign changes in the 
interior of Aq = Aq. Therefore, degQ n ,o > ^Vi,o ~ 1- If degQ ni o = iVi.o — 1 we conclude with the 
initial step. 

Suppose that degQn,o > ^i,o (including the possible case that degQ nj o = oo). It is easy to see 



that An,o(z) — A n ,o{z), so the zeros of Q n ,o come in conjugate pairs. Therefore, we can choose 
Nifl (or Nifi + 1 if necessary) zeros of Q n ,o in such a way that the monic polynomial Q* with 
this set of zeros has constant sign on A} (A} n Aj = 0). Notice that 



A n ,a 

is analytic in the indicated region and 



€ H(C\A\) 



: A "° =o(\ ) . p = V U I 



-n.O 
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In this paper the symbol £>(■) always refers to z — > oo. From JHJ, we get 

0= / x u A a ,i(x) u 7 i/ = 0,...,JV ll i-l. 

This implies that >4 n .i has at least TVi.i zeros on A}. According to Lemma [2] this linear form 
can only have N% i — 1 zeros on this interval. Consequently, our initial assumption is false and 
degQ„, = Ni :0 - 1. 

Suppose that we have proved that for some j € {0, . . . , mi — 1}, deg Q a ,j — Nij — 1, all its zeros 
are simple and lie in the interior of A j . Let us show that then, (|10p and the statement concerning 
the zeros are valid for j + 1. 

Indeed, the induction hypothesis implies that 

^ieH(C\A) +1 ), £^ = o(-L) , v = 0,...,N ltj+1 -2. 



From i®, it follows that 



da) +l {x) 



0= / x v A n , j+1 (x) J+1V / , i/ = 0,...,JV lii+1 -2. 

We have obtained (flO|) for j + 1. 

Formula (fT0|) for j + 1 implies that Qnj+i has at least iVij+i — 1 sign changes in the interior 
of Aj +1 . If degQ nj+ i = Nij+i — 1, we have finished the proof (for example, this is the case when 
j + 1 = mi because A n . mi = a„ ;mi ). Let us suppose that degQ nj -+i > Ni^+i (including the 
possible case that degQ nj+ i = oo, and of course j < m x — 2). Since A n j+i(z) = A n ,j+i(z), we 
can choose iVij+i (or 2Vij+i + 1 if necessary) zeros of Q n j+i so that the monic polynomial Q^j+i 
with this set of zeros has constant sign on Aj +2 . Then 

^^(c\A] +2 ), ^F^ ±i = °(^) ■ " " 

Using ([8]), it follows that 

>(*) 



0= x »A nJ+2 {x)-^ 



■ ^1 ,J+2 - 1 



1J + l( X )' 

This implies that -4 n .j+2 has at least N\ t j+2 zeros on Aj +2 . According to Lemma [2] this linear 
form can only have N% j+% — 1 zeros on this interval. This implies that our initial assumption is 
false; therefore, degQ nj +i = iVi^+i — 1 as stated. 

Suppose that the interval I contains two zeros x\, X2 of Q n ,j', that is, of A n ,j. According to (|10p 



•A. n ,j{v) (x-%i)(x-x 2 )daj(x) 



= , v = 0, . . . , Ni d ~ 2. 



(x - Xi)(x - X 2 ) Qn,j-l(x) 

The function A n j(x)/(x — x\){x — x%) has N\j — 3 sign changes on supp(crj) while the measure 
(x — x\){x — X2)d<jj(x)/Q n ,j-i(x) has constant sign on supp(cj). This is impossible because of 
the number of orthogonality relations. 
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Formula (fTTj) follows from (J9J) since for any q, degq < iVij+i — 1 

qAn ' j e H(C \ A) +1 ) , = o (i 



With this we conclude the proof. □ 

We need to produce additional orthogonality relations. In the second part of this section, we 
make use of some transformations employed in [llj . Let us define recursively the following functions 



(12) -An,-i-i(*) = / ^^-da](x), j = 0, 

J z X 



, m 2 . 



Proposition 3. Let n e Z™ 1+1 (.) x Z™ 2+1 («), n hmi > 1, S 1 = Af(a%, . . . , cr^J, and S 2 
N{<Tq, . . . , <t^j ) 6e given. Then, for each j = 0, . . . , 7712 



(13) J x v A a -j(x)dSj k (x) =0, k=j, ...,m 2 , ^ = 0, . . . ,n 2 ,fc - 1. 

Proof. When j = the statement reduces to the relations ii) which define ^4 n ,o- If m 2 = we are 
done. Therefore, let us assume that 7712 > 1, that (fl3|) holds for some j <E {0, . . . , m 2 — 1}, and 
prove that it is also satisfied for j + 1 . 

Fix j £ {0, . . . , m2 — 1}, k S {j + 1, . . . , 7712}, and v £ {0, . . . , 7i 2 .fe — 1}. Using the definition of 
An,— j-i, Fubini's theorem, and the induction hypothesis, we obtain 



•'■"--l,,. ; i(-' )ds» +lifc ( !B )= / a" [ ^Mda*(t)d S z +hk (x) 

• 4 - :!/) / = ^f^ dfl i+i,*( a O do i(*) = 



p„(t)A„,- J (i)do-|(t) - J t»A n ,-j{t)ds 2 hk {t) = 
since p u is a polynomial of degree < n 2l fc ~~ 2, and 7i 2 j+i > i 2 ,fe- D 

For j = 1, . . . , 7712 + 1, let Q n ,-j be the monic polynomial whose zeros are those of -An.-j in the 
region C\ A 2 ^ counting multiplicities. As we did before, in the hypothetical case that A n .~j had 
infinitely many zeros in the specified region, then Qn,-j denotes a formal infinite product. 

Taking linear combinations of the relations (|13p . we obtain 

B n2i j(x)A n ,-j{x)do-j(x) = 0, j = 0, . . . ,771 2 , 

where B n2 j is an arbitrary linear form of type 

Bn 2 ,j{x) = ^2b k (x)Sj +l k (x), degb k < n 2 ,k - 1- 
k=j 

Using Lemma [2l it follows that -4 n ,-j has at least N2J sign changes on A 2 , where 

N 2 ,j = Na t j(n) = 7i 2j H h 7i 2 , m2 , j = 1, . . . , m 2 . 

Consequently, degQ^-j > N 2 .j,j = 0, . . . , 77i 2 . Recall that for j = we proved in Proposition [2] 
that degQ nj o = — |n 2 | = |rii | — 1, that its zeros are simple, and lie in the interior of A 2 , = Aj. 
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Proposition 4. Let n e Z" il+1 («) x Z" l2+1 («), m, mi > 1, S 1 = AT(cr^ . . . , cr^J, and S 2 = 
A/"((Tq, . . . , c^ i2 ) 6e given. Then, degQ ni -j = N 2 ,j,j = 0, . . . , m,2, all its zeros are simple and lie 
in the interior of A 2 , and Q n ,-m 2 -i = 1- If I denotes the closure of a connected component of 



Aj \supp(cr|) then Q n ,-j has at most one zero in I. Moreover, 

f da 2 (x) 

(14) / x v A a ,-j{x) jy ' = 0, v = 0,...,N 2<j -l, j = 0,...,m 2 , 

and for any polynomial q, degq < N2J—1, 

(15) q(z)A a ,- j (z) = r q(x)A a ,- j+1 (x) da^jx) . ..,m 2 + l. 

Qn- 3 (z) J Q n -j(x) Z-X 2 

Proof. Fix j € {0, ... , 777.2}- From (11311 we have that for each g, deg 5 < 7J2J , 

q(z) - q(x) 



z — x 



-A n -j{x)da z Ax) = 0. 



It follows that 



A^-i-liz) = 4r / — ^(^(x) = O (1/Z" 2 - + 1 ) , Z ^ OO. 

We have shown that deg Q ni _j_i > N<2j+i(N2 >m2 +i = 0). The zeros of Q n) _j_i come in 
conjugate pair since A n ,-j-i is also symmetric with respect to the real line. If degQ nj _j_i > 
N 2 1+1 take JV2 j+i + 1 (or A r 2 j _)_ 1 +2 if necessary) zeros from Q n .-j-i so that the monic polynomial 
Qn w hh these zeros has constant sign on A 2 . If deg Q n .-j-i = N 2 ,j+i take Qn-j-i = 

Qn,-j-l- 

Therefore, 



-^n.-j-i _ Q 



and 



Qn ,-j-l 

Z -4.n.— j — 1 



(l/z n2 ^ +dcgQ -- J -i +1 ) e W(C \ A 2 ), 



Qn.-j-l 

Using ([5]), we obtain 



0(1/O eH(C\Aj), u = 0,...,n 2J + degQ* nt _ j _ 1 -l. 



da 2 (x) 

0= x A n ,-j{x)y—^ — ^ = 0, . . . , n 2 ,j + deg Q n - 1. 



This formula implies that A n .-j has at least n 2 ,j + degQ* > N 2 j sign changes on A 2 . In 
particular, we have proved that if for some j € {0, . . . , 777,2}, deg Q n ,-j-i > N 2 j+i then deg Q n ,-j > 
N 2 ,j- Going downwards on the index j we would obtain that deg Q n .o > -^2,0 = l n 2| = |ni| — 1 which 
is false according to Proposition [2] Consequently, for all j G {0, . . . , 7712}, deg Q n ,-j-i = N 2 .j+i 
(in particular, Q n .-m 2 -i = !)• Hence, Q* = Q n ,-j-i and (|14ll follows. The proof that / 

contains at most one zero of Q n ,-j is the same as in Proposition^ 

Now, fix j 6 {1, . . . ,TO2 + 1}. Notice that for any q,degq < N^ j—i, 



Qn-j Q11, 

Using ©, mi) readily follows. □ 
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3. Interlacing properties 

Fix a vector I := (hih) where < l\ < mi and < I2 < m.2. We define the multi-index 
n' := (ri! + e' 1 ;n 2 +e' 2 ) = (n^jiig ), where e li denotes the unit vector of length m.j + 1 with all 
components equal to zero except the component (Zj + 1) which equals 1. In this section it is always 
assumed that both n and n l belong to Z™ 1+1 (») x Z™ 2+1 (»). 

Fix real constants A, B such that |^4| + \B\ > and define 

G n j := AAnj + BA n i^ v j = 0, . . . ,TOi. 

Since dega^^ = degan^ + 1 it is obvious that G n j ^ 0, j < l\. In particular, this is always true 
for Q n fl. 

Lemma 4. Assume that A,BG M, |A| + \B\ > 0, and 7ii. mi > 1. Then for all j £ {0, ...,mi} 
such that m t j > 2, deg^4a n j + Ba n i j > nij — 2 and (J nj - ^ 0. 

Proof. Assume that there exists j € {0, . . . ,mi} such that riij > 2 and deg 4o nj +Ba n i j < ni t j— 3 
(rii.j — 3 = — I means that Aa nj + Ba n i j = 0). Then nj 1 — 2e 3 G Z" il+1 (*), where e J denotes 
the mi + 1 dimensional unit vector with all components equal to zero except the component j + 1 
which equals 1. According to LemmaOthc linear form G n ,o has at most |ni| — 2 zeros on Aj, but 
Gn,o satisfies the same orthogonality relations ((2|) as A n ,o and, therefore, it has at least |ni| — 1 
sign changes on this interval. This contradiction implies the statement. □ 

From this lemma it follows that if n\. mi > 2 then Q n j ^ 0, j € {0, . . . , mi}. 

Lemma 5. Assume that A, B £ K and £/ nj = AA n j + BA n i j ^ 0, for some j G {0, . . . ,mi}. 
If j < ^1 then Q n _j has at most Nij zeros, counting multiplicities, on any interval disjoint from 
Aj + i(A^ i+ i = 0). If j > l\ then Q n j has at most Nij — 1 zeros, counting multiplicities, on any 
interval disjoint from Aj +1 . 

Proof. We have 

mi 

Gn,j{z) = ^2(Aa nyk (z) + Ba n i M {z))s) +lk (z), 
k=j 

where dega„ ifc = n\,k - 1 and dega n! fe = n[\ k - 1. By LemmaH (1, sj +1>J - +1 , . . . , s] +1 m J forms 
an AT-system with respect to {n±j, ■ ■ ■ !?V mi ) on an y interval disjoint from Aj +1 , and the result 
follows immediately. □ 

Notice that for each j £ {0, . . . , mi}, Q n j is a real function when it is restricted to the real line. 

Proposition 5. Let ni :ini > 1. Assume that 4,J]el, \A\ + \B\ > 0, and let k = maxjfc' : G n ,k' ^ 
0} < mi. Then, k > l\ and Q-a,j = 0,/c < j < m\. If j < l\ then Q n .j has at most Nij zeros 
in C \ Aj + i, counting multiplicities, and at least Ni t j — 1 sign changes in the interior of Aj. If 
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l\ < j < k then G n ,j has at most Nij — 1 zeros in C \ Aj +1 and at least Nij — 2 sign changes in 
the interior of Aj . Therefore, all the zeros of Gn.j in C \ Aj +1 are real and simple. 

Proof. If j < h, then deg a n i tl > dega n ,i! and Gn,j ^ 0. Consequently, fc > l\. Obviously, from 
the definition of k, Gn,j = 0, k < j < mi. 

Assume that G n ,j,j < 'lj has at least N\j + 1 zeros in C\ Aj +1 , counting multiplicities. Select 
Ni j + 1 or Ni j + 2 zeros of Q n j which are symmetric with respect to the real axis, and let Q* ■ 
be the monic polynomial whose zeros are those prescribed. If j < l\ then 

From ©, it follows that 

/G?<7.Li (x) 
x v G^ j+ i (x) Q i +1 ^ , i/ = 0, . . . , N 1J+1 . 

These orthogonality relations imply that Gn,j+i has at least iVi.j+i + 1 zeros on Aj +1 . Since 
Gn.j+i ^ we obtain a contradiction with LemmaE] 
If j = l\ and j < k, then 

Ttf-°(?)- "- ***-»• 

Arguing as before, it follows that Gn,h+i has at least Ni^+i zeros on Aj +1 , contradicting Lemma 
[5j If j = l\ = k then (Jn^+i = and Gn.h = ^4a n ,2i + Ba n i j is a polynomial of degree at most 
"■l.ii < + 1 and thus it is identically equal to zero which is impossible. Consequently, when 
j < h, Gn.j has at most iVij zeros in C \ Aj +1 counting multiplicities. 

Let l\ < j < k and assume that G n ,j has at least N\ t j zeros in C\ Aj +1 , counting multiplicities. 
If j = mi we get immediately a contradiction because in this case Gn. mi is a polynomial of degree 
at most Ni, rni — 1. If l\ < j < m\, then there exists a polynomial Q* • with real coefficients and 
degree at least iVij such that 

^ i = °(?)- "=° 

This implies that G n ,j+i has at least Nij+i zeros on Aj +1 contradicting Lemma [SJ 

Now, let us analyze the sign changes. Notice that Gn,o ^ 0. Assume that G n ,o has N < Nx^q—I = 
| ni| — 1 sign changes on Aq = Aq, choose a nonzero linear form 

7712 

Bn 2 (z) = ^ b n 2 ,k( z )^i A z ) > deg6 n2jfc < n 2 ,fe - 1 , k = 0,...,m 2 , 

fc=0 

such that B n2 has a zero at each point where Gn,o has a sign change, and a zero of order |n.2 1 — 1 — N 
at one of the extreme points of Aq. By Lemma EJ B n2 has at most |ri2 1 — 1 zeros on Aq. Thus, 
B n2 has exactly those zeros prescribed. By definition, 

J B Yl2 {x)Gnfl{x)dal{x) = 0, 
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which contradicts the fact that B n2 (x)G n ,o(x) has constant sign on Aq. 

Let us prove by induction that for all j < lx, Q n j has at least N\,j — 1 sign changes in the 
interior of Aj. For j = this was proved above and if l\ = we are done. Let us assume that for 
some j < h, Gn,j has at least Ni j — 1 sign changes on Aj, and let us show that G n ,j+i has at least 
Nij +1 — 1 sign changes on Aj +1 . 

Let Q* ■ be a monic polynomial whose zeros are Nij — 1 points where Q n j has a sign change. 
Then 



(?). 



C'(^). i/ = 0,...,JV lji+1 -2. 



Using ([5]), this implies that 



f da],-, (x) 
= J x»G n . J+1 (x) Q i +1 ^> , v \:. s . : 2. 

Thus, Gn.j+i has at least iVij+i — 1 sign changes in the interior of Aj +1 as claimed. 

Finally, we prove that Gn,j, h < j < has at least Nij — 2 sign changes in the interior of Aj. 
Let Q n ; be a monic polynomial of degree — 1 whose zeros are points where Q-a t i t changes 
sign in the interior of A^ , then 

From here we get orthogonality conditions that imply that £/ n ,ii+i has at least N\j 1+ \ — 2 sign 
changes in the interior of Aj i+1 . One proceeds the same way until we arrive to j = k. 

From the upper bound on the number of zeros and the lower bound on the number of sign 
changes it follows that all the zeros are simple and lie on the real line. □ 

Let j 6 {0, . . . , 77i2 + 1}< Given two real constants A, B, we define 

Qn, j ■= U„ , • ISA,,, r 

Thus, by (HU, 

(16) g n .-3-i(z)= [ Gn '- j{x) da<j(x), j = 0,...,ma. 

J z — X J 

If \A\ + \B\ > then G nfi ^ and from ((TSJ) it follows that g n< _j for all j G {1, . . . , m 2 + 1}. 

Proposition 6. Let A, B € K, |A| + |S| > 0. For ewery j 6 {1, . . . , 7712}, £/n,-j Aas «i Tnosi A/2J + 1 
zeros 077 C \ A|_ 1; counting multiplicities, and at least N 2> j sign changes in the interior of A^. 
Hence, all the zeros of G n ,-j on C \ A^_j are real and simple. 

Proof. Let j S {0, . . . , 7712}. By ([13)) we know that 

a^A,! ,-j{x)ds 2 j}k {x) = 0, k = j,...,m 2 , v = 0, . . . ,n% k - 1 . 
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Since n<2,,k < n 2k> ^ follows that 

(17) J x u Q n -j(x)ds 2 jk (x) = Q , k = j,...,m 2 , v = 0, . . . , n 2 ,k - 1 ■ 

Using the same arguments employed in the previous section to show that A n ,-j has at least As- 
sign changes in the interior of A^, one obtains the same conclusion for G n .-j- 
If q is a polynomial with degq < n 2 .j, then from (|17| we have 

g(z) - g(x) 



z — x 



-Gn,-j(%)d(Tj (x) = 0. 



Hence, for every j S {0, . . . , 7712}, 

Assume that for some j G {0, . . . , m2 — 1}, G n ,-j-i has at least N 2 j + i + 2 zeros, counting multi- 
plicities, on C\ A|. Select at least N2J+1 + 2 zeros of Gn,-j-i, symmetric with respect to the real 
axis, and denote by Q* the monic polynomial whose zeros are the points selected. Then, 

ZVQ »,-i-l =0 (\), V V 2J • 1 . 



Ql,- 



As before, this implies that G n ,-j has at least iVjaj + 2 zeros in the interior of A|. Going downwards 
on the index j, we obtain that Gn,o has at least A^o + 2 = A^i o + 1 zeros, which is impossible by 
Proposition O Therefore, for all j € {1, ... , m 2 + 1}, £/n,j has at most N 2 j + 1 zeros in C \ A|_j^ 
and, therefore, they must be real and simple. □ 

Theorem 3. Let n, n' e Z™ 1+1 (») x Z™ 2+1 (•), m, mi > 2. T/ien, /or a/Z j E {-m 2 , . . .,mi} £/ie 
zeros of A n j and A n i j interlace; that is, between two consecutive zeros of A n j there is one zero 
of A n i j and viceversa. 

Proof. Since n.i, mi > 2, from Lemma |4] we know that for all j € {0, . . . , mi} and for all A, B real 
such that \A\ + \B\ > 0, the linear form Gn,j is not identically equal to zero. This is always true 
for j £ {— m 2 , . . . , —1}- Therefore, from Propositions [5] and [6] we know that for all real A,B, such 
that \A\ + \B\ > the zeros of Gn,j,j € {— ^2, ■ ■ ■ , ^i}, are real an simple. This is the basic fact 
we will use in the proof. 

Fix y G K \ Aj +1 . It cannot occur that A n ,j(y) = A n ij(y) = 0. If so, y would be a simple zero 
of A nd and A n i tj . Thus, A'^iy) f and A' alJ (y) f 0. Take A = 1 and B = -A^ d {y)/A! nlj {y) 
and consider Gn.j = AA n ,j + BA n ij. With this choice of A and B, we have 

Snjfo) = 0, 

and we obtain a contradiction because the zeros of Gn.j are simple. 

Now, taking A = A n ij(y) and B = —A n j{y), we have that \A\ + \B\ > 0. Since 

Ai< ,j (y)A n ,j (y) - A n j (y)A n i j (y) = 0, 
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and the zeros on R\ of A a i j(y)A n j(x) — A n .j(y)A Il i j(x) with respect to x are simple, it 

follows that 

But A n i j(y)A' n j(y) — A„.j(y)A' a i -(y) is a continuous real function on K \ Aj +1 in y so it must 
have constant sign on each one of the connected components of R \ Aj +1 . In particular, its sign 
on Aj is constant. 

Evaluating A n ij(y)A' n j(y) — A n ,j(y)A' nl -{y) at two consecutive zeros of A n ij, since the sign 
of A' n i . at these two points changes, the sign of A n j must also change. Using Bolzano's theorem 
we find that there must be an intermediate zero of -4 nj -. Analogously, one proves that between 
two consecutive zeros of A n j on Aj there is one of A n ij. Thus, the interlacing property has been 
proved. □ 

4. Asymptotic distribution of zeros 

Let {/i;} C A4(E) be a sequence of positive measures, where £ is a compact subset of the 
complex plane and [i G A4(E). We write 

* lim/i; = ^ , 

if for every continuous function / G C(E) 

lim J fdfii = J fdfi ; 

that is, when the sequence of measures converges to \x in the weak star topology. Given a polynomial 
qi of degree I > 1, we denote the associated normalized zero counting measure by 

q t (x)=0 

where S x is the Dirac measure with mass 1 at i (in the sum the zeros are repeated according to 
their multiplicity). 

Lemma 6. Let E C C be a compact set which is regular with respect to the Dirichlet problem and 
<f) a continuous function on E. Then there exists a unique pZ £ M.\(ET) and a constant w such that 



V»{z) + cj>{z) 



<w, z € supp(/i) , 
> w, z e E . 



If the compact set E is not regular with respect to the Dirichlet problem then the second part 
of the statement is true except on a set e such that cap(e) = 0. Theorem 1.1.3 in [27j contains a 
proof of this lemma in this context. When E is regular, it is well known that this inequality except 
for a set of capacity zero implies the inequality for all points in the set. ~p is called the equilibrium 
measure in the presence of the external field <j> on E and w is the equilibrium constant. 
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In order to determine the asymptotic zero distribution of the polynomials Q n .j we use the 
following lemma. Different versions of it appear in [5], [10] , and [28]. In [ID], it was proved 
assuming that supp(cr) is an interval on which a' > a.e. We wish to preserve this more restrictive 
condition for stronger results in section[BJ Theorem 3.3.3 in [55] and Theorem 1 in [5], do not cover 
the type of external field we consider here. So, we will sketch a proof. 

Lemma 7. Let a 6 Reg. supp(er) C R, where supp(cr) is regular with respect to the Dirichlet 
problem. Let {</>;}, I G A C Z+, be a sequence of positive continuous functions on supp(cr) such that 

(18) l ^h log wtxT\ = ^ x)> -° > 

uniformly on supp(cr). By {qi},l G A, denote a sequence of monic polynomials such that degqi = I 
and 

(19) J x k qi{x)(j)i{x)da{x) = 0, k = 0, . . . , I - 1. 
Then 

(20) *hm/ig i =/I, 
and 

1/21 

1211 lim( / K//(.c)| 2 ^(x)da(x) " 

where ~p and w are the equilibrium measure and equilibrium constant in the presence of the external 
field (f> on supp(cr). We also have that 

l/i 



(22) Urn I lw ^ 2 " ] = exp(u;- V(z)), AC C C \ Co(supp(a)). 



\oM\ 

,1/2 i, 

mm \\ej 

Proof. On account of ([T5]) and Lemma [BJ it follows that for any e > there exists Zq such that for 
all I > Iq, I G A, and z G supp(/j) C supp(cr) =: E 

7 lQ g „ ^ 4 l0 S 17771 ^ ^ + e < w - + ^ 

1 \\pi4>¥ \\e 21 

1 /2 

where € A, is any sequence of monic polynomials such that degp/ = I and ||p/<?y We = 

1 /2 

maxzeE |(pj0j' )(z)|. Hence, 

ui(z) := V^(z) + y log < w + £, zGsupp(Ti), Z > Z - 

' lb/0/ IIb 

Since it; is subharmonic in C \ supp(/l), by the continuity and maximum principles, we have 

ui(z) < w + s, z G C, I > Iq. 

In particular, 

u l (o ) = -\og———<w + e. 

1 \\Pi<t>i We 
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The last two relations imply 

(23) limsup ( I < exp (w - V^iz)), K C C, 
and 

(24) limmf Wvi^fi > ex P (-«>)• 

In particular, these relations hold for the sequence of polynomials {qi}, I S A. 

Let ti be the weighted Fekete polynomial of degree / for the weight e - ^ on supp(cr) and \a\ be 
the total variation of a. From the minimality property in the Li norm of qi, we have 

hiti /2 h ■■= (J \m(x)\ 2 Mx)da(x)^ ' < \Ml /2 \\ 2 < ^Wm^We < 

Then, using (JTSJ) and Theorem III. 1.9 in [27\, we obtain that 

(25) hmsup||^ 1/2 ||^<e-™. 

Since supp(cr) is regular with respect to the Dirichlet problem, Theorem 3.2.3 vi) in [28] yields 



hmsup l — — < 1, 



leA \\\Qi<f>ih. 

which combined with l[24|) (with pi = qi) and (|25[) implies 



(26) lim ^%^\ =1. 
Thus, we obtain dHJ) since ([21, (EHJ), and ([H give 

(27) limsup |M ; 1/2 ||^ = limsup IMpU^ = e~ w . 

leA zeA 

All the zeros of qi lie in Co(supp(cr)) C M. The unit ball in the weak star topology of measures 
is compact. Take any subsequence of indices A'cA such that 

*/™A*g ! = MA'. 

Then, 

lim y]og|<B(*)| = - lim / log— ^/^(x) = -V^'(z), JC C C \ Co(supp(a)). 

This, together with (HU) and (J22]) (applied to {qi},l G A'), implies 

(yv-v^')(z) < 0, 2 e C\Co(supp((j)). 

Since y M — V"***' i s subharmonic in C \ supp(/l) and (V^ — V /a ' a ')(cxd) = 0, from the maximum 
principle, it follows that = V^ A ' in C \ Co(supp(<r)) and thus /ia< = jl. Consequently, (f^Uj) 
holds. (HOD and (HU) imply □ 
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Using Lemma O we can obtain the asymptotic limit distribution of the zeros of the polynomials 
Qn,j, j = -m2i • ■ ■ j Tn±. At this point, let us make a slight change of notation. In the sequel, 

A 3 = A), (jj = a], j = 0,l,...,mi, 
A, = A 2 _ } , trj = o^.j, J = 0, -1, ... , -m 2 , 

and 

' 2Vij(n)-l, j =0,1..., mi, 
N2,-j{n), j = 0, -l,...,-m 2 . 
According to Propositions [2] and [H for all j — —rri2, ■ ■ ■ , mi the zeros of Q n ,j are all simple, lie 
in the interior of Aj, and total N a j points. 

Theorem 4. Let A = A( Plfi , . . . ,Pi, roi ;p 2 ,o, ■ ■ ■ ,P2,m 2 ) C Z™ 1+1 (.) x Z™ 2+1 (.), (5 1 , S 2 ) G Reg, 
S 1 = N"((Tq, . . . , c^j), and S 2 = J\f(ag, . . . , erjL) 6e given. Then 

(28) *lim//Q =J[. j = -77i2, • ■ ■ 

neA J 

where ~p — ~p(C) € Mi is the vector equilibrium measure determined by the matrix C in Q on i/ie 
system of compact sets Ej = supp(erj), j — 0, . . . , mi, = s\xpp(cr_j),j = — m 2 , ■ ■ ■ , 0. Moreover, 

\ V2|ni 



(29) ia / Era " 1 *" --l-E^ 



k=j 



where the lo£ denote the corresponding equilibrium constants. 

Proof. The unit ball in the cone of positive Borel measures is weak star compact; therefore, it is 
sufficient to show that each one of the sequences of measures {/iq n }, n € A, j = — m 2 , . . . , mi, 
has only one accumulation point which coincides with the corresponding component of the vector 
measure ~p(C). Let A' C A be a subsequence of multi-indices such that for each j = — m 2 , . . . , mi 

nGA' 

Notice that \i§ E A^i(^j), j — — 7712, . . . , mi. Therefore, 

(30) lim IQnjiz)] 1 ^ = eM-PjV^(z)), 

neA' 

uniformly on compact subsets of C \ Aj, where Pj = lim ne A' AT n .j/|ni|. 

Because of the normalization adopted on a n ,mi, An, mi = Qn,m t \ consequently, when j = mi, 
(fTUl) takes the form 

x v Q B , mi (x) J WmA ^\, =0, i/ = 0,... l JV n , TO1 -l. 
(By |cr| we denote the total variation of the measure a.) According to (|5U|) 
lim — !— log |Q», mi -i(aO| = (*) > 
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uniformly on A mi . Using Lemma [7j it follows that ^rn\ is th6 unique solution of the extremal 
problem 

P m i = w mi , x £ supp(u m , ) , 

(31) V^(x)-^^-V^- 1 (x) 1 



and 



2 ^ \ V2JV„, mi 



m lira I / lo Qn ' TOl( ^ ^ KJfo) 



n 6 A ' W |<3n,mi-l(^)| 

Let us show by induction on decreasing values of j, that for all j £ {— TO2, ■ • ■ , mi} 



(33) I ™ W _ilv«.-. W -^V»« W + %V tfl = * 



(34) Hm [ / S nA f J^'&h ^m 



where P_„ l2 _i = P mi+ i = 0, and 

l/2N nd 

,7U,lfVl 1 

n eA' y |Q nj -_i(x)| |<3n,i(a;)| 

where Q n ,-m 2 -i = 1- F° r J = m i these relations are non other than (f3"Tj) - (|3li|) and the initial 
induction step is settled. Let us assume that the statement is true for j + 1 £ {— + 1, . . . , mi} 
and let us prove it for j. 

It is easy to see that the orthogonality relations (fT0|) and (| 14[) can be expressed as 



u n M lgnJj^gMnj(g)| d\oj\{x) _ 

\Qn,j{X)\ \Q ni j-x(x)Q n<j+1 (x)\ 
On account of pip and (|15[) taking q = Q n ,j+i, this can be further transformed into 

r »Q M ( I \A nj+1 (t)\d\a J+1 \(t) \ dlajKx) = 

^ [) [J \Q aJ (t)\ \Q nJ+1 {t)\ \x-t\ ) \Qnj-l(x)Q~j+l(x)\ ' 

for v = 0, . . . , N a j - 1 . 
Relation ([3H)l implies that 

(35) hm -^loglQ^i^Q^+^x)! =-%i^- 1 (x)-%i^+ 1 (^), 

uniformly on Aj. (Since Q n .-m 2 -i = 1> when j = — we only get the second term on the right 
hand side of this limit; that is, P_ m2 _i = 0.) 
Set 

Onj+l(*) l^nj+lWI 



l2 , J+1 wi^ nj+i (:;, d|(Tj+i|W N 



IQnjWI |Qn,j + l(t)| 



It follows that for all x £ A,,- 



1 < [ Ql J+ i(t) \A n , J+1 (t)\d\a J+1 \(t) < 1 
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where < Sj +1 — inf{|a; — t\ : t S Aj+i,x S Aj} < max{|x — 1\ : t G Aj + i,x G Aj} = < oo. 
Taking into consideration these inequalities, from the induction hypothesis, we obtain that 

(36) km ( ( Q ^ d ^) 1/2NnJ e e -^+^, 

*6A' ^7 |Qnj(*)| |Qnj+l(*)| | a? — *| / 

Taking ()35[) and (|36[) into account, Lemma [7] yields that fij is the unique solution of the extremal 
problem ([55)1 and 



ton f / / 9^+1® \^j+i(t)\ d\a J+1 \(t) Ql^dW.Kx) \ 1/2jV - J = ^ 



According to (fTTj) and (fT5| with g = Q n ,j+i 

x e A 



i f Qnj+iW dfo+iKt) _ l-An,i(a;)| 



|Qn,i+l(a;)| 7 IQnJ^I |<3n,i+l(*)| \x - t\ \Q n j{x)\ 



3 ' 



which allows to reduce the previous formula to (|34[) thus concluding the induction. 

Now, we can rewrite (I33[) multiplying through by Pj and taking the constant term on the left 
to the right to obtain the system of boundary value equations 



(37) P/y^ (a:) - J " J V^-i _ 3 3+ V ^ +1 (^) 



> , x e Ej , 



for j = — 77i2, . . • , mi, where 

(38) = P/wj - P 3 P 3+1 uj 3+1 . 

The terms with P_ TO2 _i and P mi +i do not appear when j = — rri2 and j = mi, respectively. By 
Lemma [TJ (/z_ ma ,...,fi mi ) = (jl_ m2 , . . . , 7J mi ) and (w^_ TO2 , . . . , u/ mi ) = (o/ TO2 , . . . , w&J for any 
convergent subsequence showing the existence of the limits in (|28p as stated. 
Notice that implies that 

■ i -f/ i> ( ?Jo nj i1! -'i^iwV /2 ""'° 

»6A' W |Q nJ _i(a;)| \Q aJ (x)\ 



On the other hand, from (|38p it follows that P mi w mi = w^/P™, when j = mi. Suppose that 



P i+ iWj + i = Sfc 7 +i pS J : + 1 e {~ m 2 + !,•••, mi}. Then, according to 



, ,/•» mi It 



3 k=j 

and ([29]) immediately follows. □ 
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5. Proof of Theorem 1 

Here, we maintain the change of notation introduced in the previous section. Theorem [T] is a 
consequence of the following more general result. 

Theorem 5. Let A = A(pi, , . . . ,Pi >mi ;p 2 ,o, ■ ■ ■ ,P2,m 2 ) C Z™ 1+1 (.) x Z™ 2+1 (.), (S\ S 2 ) e Reg, 
S* 1 = J\f((TQ, . . . , c^jj, and 5 2 = M(cr 2 , . . . , c^ 2 ) 6e given. Let {A a j}, n G A, j = — m 2 — 1, . . . ,mi, 
6e i/ie sequences of "monic" linear forms associated with the corresponding mixed type orthogonal 
polynomials. Then, for each j = —777,2 — 1, ... , m\ 



(39) 
(A- 



lim |.4 n ,,Mr /|ni1 = Gj(z), K C C \ (A, U A j+1 ) 



m 2 — 1 



%i + l 



0), where 



(40) G,-(z) = exp ( P j+1 V^ (z) - PjV^ (z) - 2 ^ 



mi u 



Pk 



j = -m 2 - 1, . . . ,mi - 1, 



(P- 
(41) 



??l2 — 1 — ■^7711 + 1 



0) emd 



G mi (z) = exp(-P mi y^.i( z )) 



fi = ~p(C) = (jl_ m2 , . . . ,~p mi ) is the equilibrium vector measure and (w^ m2 , . . . , a;^ ) is the system 
of equilibrium constants for the vector potential problem determined by the interaction matrix C 
defined in Q on the system of compact sets Ej — supp(crj),j = 0, . . . ,mx,Ej — supp^ 2 ^), j = 
-77l 2 , ... ,0. 

Proof. If j — mi, A. n ,mi — Qn,mi an d (|28|) directly implies that 

lim|.4 n ,m l (z)| 1/|ni1 =e X p(-P mi V^(z)), KcC\A mi , 
For j G {— 77i2 — 1, ... , mi — 1}, using (fTTj) and (|15[) with q = Q n .j+i, we obtain 



Qn,j+l(z) J Qn,j(x) Qn,j + l(x) Z — X 

(Qn,-m 2 -i = !•) From (ESI)) it: follows that 

(43) lim ® nAZ ( \ =exp{P j+ iV^(z)-P j V^(z)), K C C \ (A,- U A j+1 ) 

neA Ct/n.j + ll.Zj 

(we also use that the zeros of Q n j and Q n ,j+i lie in Aj and Aj+i, respectively). It remains to 
find the |ni|-th root asymptotic behavior of the integral. 

Fix a compact set K. C C \ A J+ i. It is easy to verify that (for the definition of K 2 - +1 see proof 
of Theorem 2] above) 



Ci 



K 2 



< 



Qn,j+l( X ) Anj+iix) d(T j+ i(x) 



< 



Co 



A n,j+1 



where 



Qn,j(x) Qn,j+l(x) z-x 
minimaxi \u — x\, \v\ : z = u + iv\ : z G K,. x G A,-+i 1 
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and 



C 2 = 



i{|z — x\ : z E fC, x E Aj + i} 



< oo. 



Taking into account 



(44) 



lim 



Q a ,j+i( x ) A n , 3 +i{x) da j+ i(x) 



Qn,j(x) Q„ J + l(x) Z-X 

From (14"2"1) - (|4"4"1) . we obtain (|3"9")) and we are done. 



i/KI 



exp ( -2 £ </P* 
k=j+i 



□ 



Remark 1. Taking into consideration that the polynomials Q n .j (see Propositions [2] and [4]) and 
the functions 

Ql,j( x ) A nJ {x) dcrjix) 



)n,j-l{x) Qn,j{x) Z-X 1 

may have at most one zero in each of the connected components of Aj \ Ej , in place of ((39)) one 
can prove convergence in capacity on each compact subset K, C C \ {Ej U Ej + i). More precisely, 
for any such compact set JC and each e > 

limcap(ze/C: |-4 nj (z)| 1/|ni - G 3 \z) > s\ = 0. 

□ 



Set 



U*(z) = PjV* (z) - P J+1 V»w (z) + 2 



mi jt 



k=j+l 



Pk 



-m% — 1, . . . , mi — 1, 



and 



c/ mi (z) = p roi y^(2). 



Hence, Gj(z) = exp(-U?(z)),j = -m 2 - 1, 



We have that for j 



-TO 2 , 



nil • 

mi(P- m2 -l = Pmi+1 = 0) 



v*j+i (z) + pfv^ (z) - / ' ,/ ; ' i '• (z) 



■ w7 



From the equilibrium property (see Lemma [T] and (|37|) ). it follows that 

= 0, x E suppQ^) , 



> 0, x E E; 



Define 

Pl,j, j = 0,...,TOl, 

It is easy to verify that for j — —m%, . . . , m\ 



(45) 



U$[z) - U»{z) = 0{{ Pj -ft-i)log l/\z\), z - oo. 



In particular, Uj{z) — C/j 1 _ 1 (z) = 0(1), z — > oo, whenever p,,- = Pj-i- By assumption, p., — Pj-i < 
0,j = -m 2 ,...,wi except for p -p_ x = Pi, +P2,o > 0. 
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For all j, the function Uj — Uj_ 1 is subharmonic in C \ supp^). If pj > Pj-i, then it is 
subharmonic in all C\ supp(/Ij)- According to what was said above, when j — or pj — Pj-i, from 
the equilibrium condition and the maximum principle, we have that 17? — U?_ x = on supp(<T :) ) = 
Ej and Uj < on C \ supp(crj). In particular, in this case we have that supp(/I J ) = supp(crj). 

When Pj-i > Pj, (|15|) implies that in a neighborhood of z = oo,U? > U^_ 1 . Let jj = {z £ 
C : Uj(z) — t r j J _ 1 (z) = 0}. The equilibrium condition entails that jj D supp(7tj) and the initial 
remark indicates that jj is bounded. Consider any bounded component of the complement of jj. 
On it, Uj 1 — Uj t _ 1 is subharmonic and on its boundary — ^j-i = 0- Thus, on any bounded 
component of the complement of jj we have that Ujf < U^_ x . From the initial remark it follows 
that on the unbounded component of the complement of jj, Uj > U?_^. 

Fix j g {0, . . . , mi}. For each k £ {j, . . . , mi} define 

D{ :={zeC\ U?L l -Aj : Uf{z) < Uf(z), i = j, . . . , m u i ± k}, := C \ A roi . 



Let 



Cj(z) = vam{U?(z) :k=j,..., mi}. 



Corollary 1. Let A = A(p 1|0 , . . . ,Pi, mi ;P2,o, ■ ■ ■ ,P2,m 2 ) C Z™ 1+1 (.) x Z™ 2+1 (.), (S 1 , S 2 ) e Reg, 
S 1 = N{ul, . . . , cr} ni ), S 2 = Af(po, . . . , of ra2 ) fee groen. Let (a n ,o, a n ,i, ■ • ■ , s^mji n e A, fee ifte as- 
sociated sequence of "monic" mixed type multiple orthogonal polynomials. Then, for j = 0, . . . , mi 

(46) lira \a n ^)\ 1/M = exp(-G(z)), £ C Uj^Dj, 
and 

(47) limsup|a nj -(z)| 1 /l ni l < exp(-&-(*)), KcC\U™A. 

neA 

In particular, if Pi,o = ■ ■ ■ = Pi. mi = l/(mi + 1), te^en 

(48) lim |a„.,(z)| 1/|ni1 = exp(-C/^( Z )), KcC\ U^.A fe . 

nGA J 

Proof. For j = mi, A n , mi = Gn.rou-DmJ = C \ A mi and £ mi = J7£ . Therefore, ([46]) reduces to 
(|39|) and implies 1|47|1 . Let us prove these relations for j = 0, . . . , m\ — 1. 

The *4 n j are expressed in terms of the a n ,fc, k = j, . . . , mi, through a linear triangular scheme 
of equations with function coefficients which do not depend on n. Using this system, we can solve 
for a n j, in terms of A n .k, k = j, . . . , m±. 

Given j G {1, . . . , mi} and < i < j, we have 

r _ lV -i/ l iw \ __ f f dal{x l )---dcj}{x j ) 

(I) { a i ,...,a j){ z, J - J {z _ Xi){x ^ _ Xi) ^ ^ _ ^ , 

where (-){z) denotes the Cauchy transform of the indicated measure, and 
i 10 s [ [ da}{x l )---da){x j ) 



3 ' ' 



(iEj+l - Xi) ■ ■ ■ (Xj - Xj-i)(z - Xj) 
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Consequently, 

— {xj - Xi)da}(xi) ■ ■ ■ daUxj) 



(z - Xj)(x,+1 - Xi) ■ ■ ■ (Xj - Xj-l)(z - Xj) 

Since Xj — Xi = Xj — Xj—i + Xj—i — ■ ■ ■ — 2^+1 + x,_|_i — x,-, substituting this in the previous formula, 
we obtain 

3-1 

(49) {a),..., a}){z) = £>l) fc ->i, 4?(^h + (-!)''-><> • • • , 

k—i 

(This formula is applicable to any Nikishin system. We will use it on S 2 in the last section.) 
Using formula (I49p it is easy to deduce that (the sum is empty when j = mi) 

mi 

o^W = A^j{z) + (- 1 ) k ^(^l ■ ■ ■^}+i)(z)A nM (z). 
k=j+l 

Taking (I39| into consideration, on Dl the term containing A n .k dominates the sum (notice that 
(o\ : . . . , crj +1 ) (z)/0,zeC \ Afe) and (|46|) immediately follows. On the complement of Uj™^ ■£)][. 
there is no dominating term and all we can conclude from the previous equality is (|47|) . 

Let pi.o = • • • = Pi, (Tii = l/( m i + !)■ bi this case, on C \ U]™^Afc we have that U{f ll (z) < 
Umx-ii 2 ) < ■ ■ < Uf(z) and (gHJ) follows from (46]). □ 

6. Ratio asymptotic 

Here, we study the convergence of the sequences {Qn l ,j/Qn,j}, n S A C Z™ 1+1 (») x Z" l2+1 (») 
and of the ratio of the corresponding linear forms. We maintain the notation introduced in Section 
[H namely 

Aj = A], Oj = o-j, j = 0, 1, . . . , mi, 
A j = A 2 _ j , a j = a 2 _ j , j = 0, -1, . . . , -m 2 , 



and 

(50) jV n j = 



Set 



iVij(n)-l, j = 0,1..., mi, 
^2,-i(n), j = 0,-1, ...,-m 2 . 

Qn,j+1-An,j . -. 

n nj = — , j = —mi — l,...,mi, 

^ n >3 

(Qn,-ra 2 -i = Q n ,mi+i = 1 an d W n , mi = 1). With these notations, relations (fT0|) . (fl4| . (fTTjl . and 
(I15|) (replacing general q by Q n ,j+i and shifting the index j by —1) can be rewritten as follows 

(51) f x"Q nd (x) . l^^M^-Kx) = (/ „ v i. i = _ TO2) ..., mi! 

7 J / |Q n , J -_i(x)Q nJ+ i(x)| 

and 

( 62 , ^-.i^/ t':' '"fy",v * 

J Z-X (y n> j_i(x)Qn,i+l(x) 
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Since on the interval Aj the measure Uj and the functions Hnj, Qnj-iQn.j+i, preserve a constant 
sign, we can take their absolute values in (|5ip without altering the orthogonality relations. 
For each j = — m 2 , . . . , mi, define 

-1/2 



(53) 
Take 



K _ ( [ n i <„\ l^n,j(a:)| d\a 3 \{x) 



K a , 3 . 

An.mx+l = 1 , Kn,j = 17 , J = -m 2 , ...,171! . 



Define 

(54) q n<j = K nj (3 nj , h n ,j-i = K^M n ,j-i , 



and 

(55) dp nJ (x) 



h n _j(x)d(jj(x) 



Qn,j-l{x)Q a ,j+l(x) 

From (|5ip and the notation introduced above, we obtain 



(56) / x v Q n j(x)d\p nd \(x) = 0, ^ = 0, ...,N n j - 1, j = -m 2 , • . . , m Y , 

J Aj 

and g n j is orthonormal with respect to the varying measure |/5 nj |. On the other hand, using 
it follows that 

(57) h n ,j-i(z) = e nJ / — rf|p nj -|(x), j = -m 2 , . . . ,mi, 

where e nj denotes the sign of the varying measure p n j- 

The proof of Theorem [S] below has three steps. First, we show that for each j G {— . . . , mi} 
the sequence {Q n i j/Q n ,j} is uniformly bounded on each compact subset contained in C\ supp(tTj) 
(for all sufficiently large |ni|). Taking a subsequence of multi- indices such that all the sequences 
of ratios of polynomials have limit, we show that the limit functions must satisfy a system of 
boundary value problems. This system happens to have a unique solution from which we derive 
that all convergent subsequences have the same limit. Finally, we show that the limit functions can 
be expressed in terms of the branches of certain conformal representations of a related compact 
Riemann surface onto the extended complex plane. 

In this section, we assume that supp(crj) = Aj U ej,j — —m 2 , ■ ■ ■ ,mi, where Aj = [cij,bj] is a 
bounded interval of the real line, \oj\ > a.e. on Aj, and ej is a set without accumulation points 
in R\ Aj. We denote this writing S 1 = Af'(<r%, . . . , ct^J, S 2 = Af'(o%, . . . , of„ 2 ). 

In order to fulfill the first step, Theorem [3] would be sufficient if Aj — Aj, j = —m 2 , . . . , mi. In 
order to allow the compact sets to enter the connected components of Aj \ supp(erj), we need to 
show that the zeros falling in the intervals / (see Propositions [5] and 2]) are attracted to points in 
supp (<7j) \ Aj. In our aid comes the next result. 
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Lemma 8. Let S 1 = N'{al, . . . , cr^ 1 ), S 2 = Af'(o-%, ...,a^ 2 ) be given, and let A C Z™ 1+1 (») 
Z™ 2+1 (») be an infinite sequence of distinct multi-indices such that 

(58) sup((m 2 + l)n 2 ,o - |n 2 |) < oo , sup ((mi + l)ni >0 - |ni|) < oo . 

neA neA 



For any continuous function f on supp(<Zj) 
(59) lim / f(x)ql j (x)d\p nJ \(x) = f / f(x 



neA 



1 f bj .. , dx 



TT J aj y/(bj -x)(x-aj) 

where Aj — [aj, bj], —mi <j<m\. In 'particular, 

(60) lim £ n ,jftn,3-l( z ) = 7 ? \, =^ > /C C C \ SUpp(cr J ) , 

neA y/iz-bj^z-aj) 

where \J {z — bj)(z — eij) > «/z > 6j. Consequently, for j — —mi, . . . , mi, each point o/supp(cjj)\ 
A-,- is a Zirazi o/ zeros 0/ {Q n ,j}, nGA. 

Proof. We will prove this by induction on j. For j — mj, using Corollary 3 in [5] and the second 
condition in (I5"8"l). it follows that 



hm / /(a;)< mi (a:)|^ 737 = -/ /W 



nGA 



where / is continuous on supp(cr mi ). Take /(x) = (z — x) _1 where z € C \ supp(cr mi ). According 
to (|57p and the previous limit one obtains that 

lim En, roi /in, mi -i(2:) = 7= ===== =: ftmi («) , 

neA V( z ~ &mj(-2 - ami) 



pointwise on C \ supp(o" mi ). Since 



yn.mi (z) rf|cr mi |(x) 



- HIP ~1 vT> z e K, cC\supp(a mi ) 

d{IC, supp(cr m J) 



Z % I Qn.mi — 1 (■&) 

where d(/C, supp(<7 mi )) denotes the distance between the two compact sets, the sequence {/i n ,mi-i}, 
n 6 A, is uniformly bounded on compact subsets of C \ supp(cr mi ) and (|60[) follows for j = mi- 

Let £ S supp((T mi ) \ A mi . Take r > sufficiently small so that the circle C r = {z : \z — £| = r} 
surrounds no other point of supp(cr mi ) \ A mi and contains no zero of g n ,mi, nGA. From (|60p for 
j = m\ 

(61) Um-L [ ^-fj^-L / £^ = 0. 

neA 27Ti J Cr e n ,mi/J„, mi -i(z) Z7tz J Cr h mi (z) 

Since C is a mass point of <r mi , formula (|57p indicates that either ft. n ,mi-i has a simple pole 
at £ or Q n ,mi(C) = 0- I 11 anv case, from (|6ip and the argument principle, it follows that for all 
sufficiently large |n|,n G A, Q n . mi must have a simple zero inside C r . The parameter r can be 
taken arbitrarily small; therefore, the last statement of the lemma readily follows and the basis of 
induction is fulfilled. 
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Let us assume that the lemma is satisfied for j G {k + 1, . . . , mi}, — 7712 < k < mi — 1, and let 
us prove that it is also true for k. From (f60|) applied to j = k + 1, we have that 



n£A 



i/| (x- b k +i)(x- a k +i)\ ' 
uniformly on A& C C\supp(ofc + i). It follows that {|/i n .fc|d|cfc|}, n G A, is a sequence of Denisov type 
measures according to Definition 3 in |2_. Additionally, ({|/i n ,fcM|ffe|}, {|<3n,fc-i<3n,fc+i|}, I), n € A, 
is strongly admissible as in Definition 2 of for each I G Z (see paragraph just after both 
definitions in the referred paper). Therefore, we can apply Corollary 3 in 2J of which (|59|) for 
j = k is a particular case. In the proof of Corollary 3 of [2] (see also Theorem 9 in [4]) it is required 
that the inequality deg(<5 nj _iQ nj +i) — 2deg(Q nj ) < C holds for every neA, where C > is a 
constant. It is straightforward to check that this condition is satisfied under the assumption (|58p . 

Now we return to the induction argument. From (f59"| for j = k, and the rest of the 
statements of the lemma immediately follow just as in the case when j = mi- With this we 
conclude the proof. □ 

Now, we are ready to prove normality. 

Lemma 9. Let S 1 = Af'(<T^, cr^J, S 2 = J\f'(a 2 , a 2 mQ ) he given, and let A C Z™ 1+1 (.) x 
Z™ 2+1 (») be an infinite sequence of distinct multi-indices such that 



sup((m 2 + 1)^2,0 - |n 2 |) < 00 , 

n£A 



sup((mi + l)ni,o — |ni|) < 00 . 

neA 



Let us assume that there exists I = (Zi ; Z2) , < l\ < mi,0 < I2 < rn-z, such that for all n G A 
we have that n' G Z™ 1+1 (») x Z™ 2+1 (»). Then, for each j = — n%2, ■ ■ ■ ,T7ii, and each compact set 
K, C C \ supp(<T :) ) there exist positive constants Cj t iQC),Cj t 2()C) such that 



C jtl (fC) < inf 



zeK. 



/(*) 



Qn,j(z) 



< sup 

zeK. 



Qn'.j(z) 



Qn,j{z) 



< cm, 



for all sufficiently large |ni|,n G A. 



Proof. The uniform bound from above and below on each fixed compact subset K, C C \ Aj (for 
all n G A) is a direct consequence of the interlacing property of the zeros of Q n i j and Q n ,j- I n 
fact, comparing distances to z G K, of consecutive interlacing zeros, it is easy to verify that 



iin K!|} 



< inf 

2 J z£)C 



<9n,j(z) 



< sup 

zeic 



Qn,j{z) 



< 



max{ii2, d 2 ,} 



where di denotes the diameter of K, U A^- and d\ denotes the distance between IC and Aj . So, for 
such compact sets the assertion has been proved. 

The additional restrictions made in the lemma guarantee that the zeros of the polynomials 
Q n i j and Q n ,j lying in Aj \ supp(<jj) converge to the mass points as results from Lemma [8] Let 
IC C C \ supp(crj) and suppose that K, <~) Aj ^ 0. Notice that K. can intersect at most a finite 
number of open intervals . . ,Im forming the connected components of Aj \ supp(<7j). The 



MIXED TYPE MULTIPLE ORTHOGONAL POLYNOMIALS 31 

polynomials Q n 'j and Qn,j can have at most one zero in each of those intervals. Consequently, 
for all |ni|,n G A, sufficiently large, the zeros of Q n ij and Qnj lie at a positive distance e from 
JC. Now, it is easy to show that for all sufficiently large |ni| 



■{•■=} 



< inf 

2 J zdK 



?n,j(z) 



< sup 



Q n iAz) 



Qn,j{z) 



max{(i2, d^} 



This concludes the proof. □ 
From Lemma [S]we know that the sequences 

{ ( 9n'j/Qn,j} n6A , 3 = ~™>2, ...,mi, 

are uniformly bounded on each compact subset of C \ swpp(crj) for all sufficiently large |ni|. By 
Montel's theorem, there exists a subsequence of multi-indices A'c A and a collection of functions 
Fj , holomorphic in C \ supp(crj), such that 



(62) lim Q nfl =FP{z\ £ C C \suppfo). j = -m 2) . . . , mi. 



3 

In principle, the functions F^p may depend on A'. We shall sec that this is not the case and, 
therefore, the limit in (f6"2")l holds for n e A. First, let us obtain some general information on the 
functions F^p . 

The points in supp(<7j) \ Aj are isolated singularities of Fj . Let £ G supp(cj) \ Aj. By 
Lemma EJ £ is a limit of zeros of Q n .j and Q n ij as |ni| ~> oo,n € A, and in a sufficiently small 
neighborhood of £, for large |ni|,n S A, there can be at most one zero of these polynomials (so 
there is exactly one, for all sufficiently large |ni|). Let lim ne A £ n = £ where Q n ,j((n) — 0. From 
formula ([62]) 

hm (iz^O^g) = (z _ j?ro (js) c c (c \ ( }) □ {c} 

neA' Qn, 3 (z) 

and (z — C)fP (z) is analytic in a neighborhood of £. Hence C is not an essential singularity of F^ . 
Taking into consideration that Q n ij, neA, also has a sequence of zeros converging to £, from the 
argument principle it follows that £ is a removable singularity of which is not a zero. By Lemma 
[9] we also know that the sequence of functions \Q B ij/Qa,j\, n S A, is uniformly bounded from below 
by a positive constant for all sufficiently large |ni|. Therefore, in C\ supp(crj) the function Fj is 
also different from zero. According to the definitions of Q n ,j, Q n l j, and Propositions [2] and [4] (see 
also ([50]) ). when —l 2 < j < h, we have that degQ n ij — N n i t j = N n j + 1 = degQ nj + 1 whereas, 
for j S {—m2, ■ ■ ■ , — fa — 1} U {1% + 1, . . . , mi}, we obtain that deg Q n i j = N a i j — N n j — degQ nj . 
Consequently, when — fa < j < l\, the function F^p has a simple pole at infinity and (fP)'(oo) = 1, 
whereas, for j G {— m,2, ■ ■ ■ , — fa — 1} U {fa + 1, . . . , mi}, it is analytic at infinity and i^(oo) = 1. 

Now, let us prove that the functions Fj ,j = — mi, . . . , mi, satisfy a system of boundary value 
problems. 
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Lemma 10. Let S 1 = A/ 7 ((To, . . . , cr^J, S 2 = TV' (erg, ...,<r^ 2 ) be given, and let A C Z" ll+1 (« 

+ 



Z^ 2+ (•) 6e an infinite sequence of distinct multi-indices such that 



sup((m 2 + l)ri2,o - l n 2|) < oo, sup ((mi + l)n lj0 - |ni|) < oo. 

nEA ngA 

Let us assume that there exists I — (£15/2)1 < l\ < mi,0 < Z2 < ?n 2 , smc/j that for all n £ A we 
/iaue f/iaf n' G Z™ 1+1 (») x Z™ 2+1 (»). Then, there exists a normalization Fj , j = — m 2 , . . . , mi, 6?/ 
positive constants, of the functions Fj given in which verifies the system of boundary value 

problems 

1) ff, effCCXA,-), 

2) (Ff)'(oo)>0, j e {-/ 3 ,.--,ii}, 

(63) 2') if>(oo)>0, je{-m 2) ...,-Z2-l}U{Zi + l,...,mi}, 

3) |if } (z)| 2 7 7r -^ = 1, x e Aj , 

1(^1^1)^)1 

^ere = = 1. 

Proof. The assertions 1), 2), and 2') were proved above for the functions Fj l \ Consequently, they 
are satisfied by any normalization of these functions by means of positive constants. 
From (|5o]) applied to n and n z , for each j = — m 2 , . . . ,nii, we have 

x"Q nt j(x)d\p n>j \(x) =0, v = 0, .. . ,iV n j - 1 , 



and 



where 



x v Q n ij{x)g n j(x)d\p n ,j\(x) = 0, ^ = 0, . . . , N n i tj - 1 



? / s _ IQn,j-i(^)Qn, 3 +i(a:)| Ifen'jWI ^ / s _ ^n, j {x)da J (x) 

nJ \Qn',j-l(x)Qn',j + l(x)\ \hn,j(x)\ ' " J Q n ,j-l(x)Q n ,j+l(x) 



From (HOD and dB2J) 



(64) ]hz i g n j(x) = \(F® 1 F$ 1 )(x)\- 1 
uniformly on A.,-. 

Fix j G {— 77i 2 , . . . , — ?2 — 1} U {/1 + 1, . . . , mi}. As mentioned above, for this selection of j we 
have that degQ n i j = degQ n j = N n ,j- Due to (|64|) and |62|) . from Theorems 1 and 2 of [2], it 
follows that 

(65) lim = ItH = ^ = ^\z) , K C C\ suppfe) , 
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where Sj is the Szego function on C\ Aj with respect to \Fj_ 1 (x)Fj_^ 1 (x)\~ 1 , x G Aj. The function 
Sj is uniquely determined by 

1) SjA/SjeHiCXAj), 

(66) 2) Sj(oo)>0, 

Now, fix j £ {—h, • • ■ , h}- In this situation degQ n ; j = degQ nj + 1 = N n j + 1. Let Qnji x ) be 
the monic polynomial of degree -/V nj orthogonal with respect to the varying measure g a j d\p n j\. 
Using the same arguments as above, we have 

(67) i^^lM M± s j{zh KcC\supp(a 3 ). 

On the other hand, since deg Q n i j = deg Q* ■ + 1 and both of these polynomials are orthogonal 
with respect to the same varying weight, by Theorem 1 of 2\ and ([60]) , it follows that 

(68) hm -^-^ = — BW (,j, K C C \ Buppfo) , 

where i/?j denotes the conformal representation of C \ Aj onto {w : \w\ > 1} such that ipj (oo) = oo 
and ip'j(oo) > 0. The function ipj is uniquely determined by 

1) W) lM€ff(C\A j ), 

(69) 2) ^(oo)>0, 

3) |^(x)| = l, x£A r 

From (J62J), (67]), and (J6SJ) , we obtain 

(70) hm 9fM = (s j <p j )(z)=FP(z), /CcC\supp( ( r,). 

n£A' Qn,j{Z) 



Thus, 
(71) F® 



S 3<?o, 3 £ {-h,---,h} , 

S 3 , j e {-m 2 , . . . , -Za - 1} U {Zi + 1, . . . , mi} , 



and from (|6l)|) and (|71 [) it follows that 

(72) \Fj l \x)\ 2 1 = — , zeAj, j = -m a ,...,mi 



where 



(73) 



(S^M, je{-Z 2 ,...,Zi}, 



J 3 Yji 
3 1 



5 2 (oo), je {-m 2 ,...,-Z2-l}U{Z 1 + l,...,mi}. 
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Now, let us show that there exist positive constants Cj , j — — 7712 , ■ ■ ■ , rn\ , such that the functions 

3 j 



= CjF^ satisfy ((63)) . In fact, according to (l72j) for any such constants Cj we have that 



1 c 2 

\F' ?\x)\ 2 - — m — m = , x€Aj, j = -m 2 ,...,mi, 

where c_ m2 _i = c mi +i = 1. The problem reduces to finding appropriate constants Cj such that 

(74) J - = 1, j = -m 2) ...,mi. 

Taking logarithm, we obtain the linear system of equations 

(75) 21ogCj ■- logCj_i - logCj +1 = logujj , j = -m 2 , . . . , m 1 

(c-m 2 -i — c mi +i = 1) on the unknowns logo, . This system has a unique solution with which we 
conclude the proof. □ 

Consider the (mi + + 2)-sheeted Riemann surface 



n = [j n k , 

k— — ni2 — 1 

formed by the consecutively "glued" sheets 

H-m 2 -\ ■= C\ A_ m2 , TZ k := C \ (A fc U A fe+ i), fc= -m 2 ,...,mi - 1, 7^ mi := C \ A mi , 

where the upper and lower banks of the slits of two neighboring sheets are identified. Fix I = 
(h,k),0 < h < TOi,0 < h < m 2 . Let be a singled valued function defined on TZ onto the 
extended complex plane satisfying 

^(0(^ = ^ + 0(1), z^oof-^- 1 ), 

z Z A 

i/>®(z) =C 2 z + 0(l), z^oo^, 

where C\ and C2 are nonzero constants. Since the genus of TZ is zero, ip^ exists and is uniquely 
determined up to a multiplicative constant. Consider the branches of tp"', corresponding to the 
different sheets k = —mi — 1, ... , m\ of TZ 

^(0 :={^ ) }^_ ma _ 1 . 

We normalize so that 

mi 

(76) [] 1^(00)1 = 1, C ie K\{0}. 

fc— — mi — 1 

Certainly, there are two verifying this normalization. Since the product of all the branches 
rifeS— m 2 -i ^ s a sm gl e valued analytic function in C without singularities, by Liouville's theorem 
it is constant and because of the normalization introduced above this constant is either 1 or —1. 
If ipV) is such that C\ G M \ {0}, then 



^«(z)=^W(z), zeTZ. 
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In fact, let <f>(z) := i/jW(~z). <f> and ip^ have the same divisor; consequently, there exists a constant 
C such that (f> = Cip( l \ Comparing the leading coefficients of the Laurent expansion of these 
functions at oo'- - ' 2 " 1 "' , we conclude that C = 1 since C\ G M\ {0}. 

In terms of the branches of the symmetry formula above means that for each k = —mi — 
l,...,m 1 , 

4° :l\(A t UA w )^I 

(A_„ l2 _i — A,„ 1+ i — 0); therefore, the coefficients (in particular, the leading one) of the Laurent 
expansion at oo of these branches are real numbers, and 



(77) ^ ) (x ± )=^\x T )=^l 1 (x ± ), xeA k+1 . 

Given an arbitrary function F{z) which has in a neighborhood of infinity a Laurent expansion 
of the form F(z) = Cz k + ©(z^ 1 ), C ^ 0, and k G Z, we denote 

F :=F/C. 

C is called the leading coefficient of F. When C G R, sg(F(oo)) will represent the sign of C. 
We are ready to state and prove one of the main results of this section. 

Theorem 6. Let S 1 = A/"'(ctq, . . . , a^ ni ), S 2 = jV'(<t^, . . . , a^ 2 ) be given, and let n G A C 
Z" 11+1 (») x Z" l2+1 (») be an infinite sequence of distinct multi-indices such that 

sup((m 2 + 1)^2,0 - |n2|) < oo, sup((mi + l)ni,o - |ni|) < oo. 

n£A nSA 

Let us assume that there exists I — (£15/2)1 < l\ < mi,0 < 1% < m2, such that for all n G A we 
have that n l G Z™ 1+1 (») x Z™ 2+1 (»). Let {Qn,j}™J;- m2 , n G A, be the corresponding sequences of 
polynomials defined in section^ Then, for each fixed j G {— mi, . . . , mi}, we have 

(78) lim = *f (*)> z£ K C C \ SU PP^) 

n£A Q n! j(Z) J 

where the functions satisfying JTjjj) are 



(mi \ mi 

Proof. Since the families of functions 

{Qn',j/Qn, 3 } neA , j = -TO 2 , ...,1711, 

are uniformly bounded on each compact subset £cC\ supp(<7j) for all sufficiently large |ni|, n G 
A, uniform convergence on compact subsets of the indicated region follows from proving that 
any convergent subsequence has the same limit. According to Lemma [10] the limit functions, 
appropriately normalized, of a convergent subsequence satisfy the same system of boundary value 
problems (l63|) . According to Lemma 4.2 in [1] this system has a unique solution. 
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It remains to show that the functions defined in (|79|) satisfy (|63|) . When multiplying two 
consecutive branches, the singularities on the common slit cancel out by the Schwarz reflection 
principle; therefore, 1) takes place since only the singularities of ijrp on Aj remain. From the 
definition of ip^ it also follows that for j = —l 2 , . ■ ■ , 1%, ff has at infinity a simple pole, whereas 
it is regular and different from zero at infinity when j 6 {— TO2, • • • , — h — 1} U {li + 1, . . . , mi}. 
The factor sign in front of ([79]) guarantees the positivity claimed in 2) and 2'). 

In order to verify 3), notice that ff /ff i = sg(^f ^oo))/^.: x . Therefore, if j — —m 2 + 
1, . . . ,mi, 

|ff 0r)| a \i>f{x)\ 



on account of (f77| . For j = — 7722, from the definition and ((771 

mi mi 

- u/,W 



= 1, i e 



Am 2 

|-^ 7 '-m2+l( :Z; )l - I. — n> 2 i 

S) 



since Ilfe=— ma— l ^fc * s constantly equal to 1 or —1 on all C. The proof is complete. □ 

The following corollary complements Theorem [6l The proof is similar to that of Corollary 4.1 
in p. 



• x 



j/m 2 + l 



Corollary 2. Let S 1 = N'(pi, . . . , o^J, S 2 = Af'(a 2 , . . . ,o^ 3 ) be given, and let A C Z™ 1+1 ( 
(•) be an infinite sequence of distinct multi-indices such that 

sup((m 2 + l)n 2 ,o - |n 2 |) < oo, sup((mi + l)ni )0 - |ni|) < oo. 

n£A nSA 

Let us assume that there exists I = (Zi ; Z2) , < l\ < mi,0 < l 2 < ?7i2, suc/i i/iai for all n G A 
we ftaue t/iaf n' € Z™ 1+1 («) x Z™ 2+1 («). Lei {g„j = K B j Q n j}J^_ m2 ,n E A, &e the system of 
orthonormal polynomials as defined in \54-\ and {Kn,j}^l_ m2 ,n £ A, the values given by ([5 
Then, for each fixed j = —m 2 , ■ ■ ■ , mi, we have 

(80) lim ^ = K (0 



(81) lim |^ =K (0... K W 
and 

(82) lim = ^ K C C \ SU PP(^) ' 

n£A q n j{Z) 



whe 



re 



(83 ) K W = , c (i) 



(Ff)'(oo), je {-h,...,h}, 

^f(oo), 

( c -m,-i = c mi+i = 1) anc ^ ^ e functions ff are defined by ( |7ffp . 
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Proof. By Theorem [6l we have limit in (|64| along the whole sequence A. Reasoning as in the 
deduction of formulas (|65p and (|70[) . but now in connection with orthonormal polynomials (see 
Theorems 1 and 2 of [2]), it follows that 



lim 



q a i,j{z) J [Sj<pj)(z), j e {-h,-..,h}, 



nEA ^W 1 5,(2), iG{-m 2 ,... 1 -J 2 -l}U{Zi + l,...,mi}, 

uniformly on compact subsets of C \ supp(<7j), where Sj is defined in (|66[) . This formula, divided 
by (|65|) or (f70|) according to the value of j gives 



urn — = ,M = 



\/ c j-l c j+l 

where u>j is defined in (|73p. and the Cj are the normalizing constants found in Lemma 1 101 solving 
the linear system of equations ()75() which ensure that 



^ = Cj*j , J — ~m 2 , ■ ■ ■ , mi , 

with Ff ] satisfying |JB3J and thus given by IJ73J. Since (Fj !) )'(oo) = l,j € {-Z 2 , • • • , Zi}, and 
(i^. )(oo) = 1, j 6 {— 7712,...,— '2 — 1} U {h + l,...,mi} formula ([80]) immediately follows with 
^° as in (PI). 

From the definition of K nj - , we have that 

Taking the ratio of these constants for the multi-indices n and n' and using (|Bu| . we get (f8"Tj) . 
Formula 1(52]) is an immediate consequence of ([50]) and (|75|) . □ 

Let lcm(a, b) denote the least common multiple of two integers a and b, and define m := 
1citi(toi + l,m2 + 1), d\ := m/(mi + 1), c?2 := m/(m 2 + 1). Within the class of pairs I — (Zi ; Za) 
with < Zi < mi, < Z 2 < m 2 , we distinguish the subclass 

L := Z 2 ) ■ h = r mod (mi + 1), Z 2 = r mod (m 2 + 1) for some < r < m — 1} . 

It is easy to check that for different r, < r < m — 1, the pairs (Zi,Z 2 ) in L are different. Let 
p := (pi; p 2 ), where pi = (di, . . . , d±) and p 2 = (cZ 2 , . . . , d 2 ) have mi + 1 and m 2 + 1 components, 
respectively. By n + p we denote the multi-index (ni + pi; n 2 + p 2 ). 



Corollary 3. Let S 1 = A/Vo, cr^J, S 2 = 7V'(cr 2 , . . . , cr^J 6e given, and let A C Z" ll+1 ( 



• 1 x 



(•) &e an infinite sequence of distinct multi-indices such that 

sup((m 2 + 1)^2,0 - l n 2|) < 00, sup((mi + l)ni )0 - |rii|) < 00. 

n£A ngA 

Then, for each fixed j S {— m 2 , . . . , mi}, we have 

(84) Um%±MM = n^ ) (^ /CcC\su PP (<7,), 

n£A L^nj'W ;gL 
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(85) 

and 
(86) 



lim 



'n+PJ 



-n«s°. 



Proof. Given n 6 A and < r < m, let n(r) := n + q(r) where q(r) = (qi(r); q2(r)) is the 
multi-index satisfying 

q, 4 (r) = (A; + 1, . . . , k + 1, k, . . . , k) , r — k(m t + 1) + s , < s < to, . 



Hence, n(0) = n, n(m) = n + p and n(r) € Z™ 1+1 (») x Z™ 2+1 (») for every r. 
We have 

Qn+p,j(^) _ Qn(r+1), 
Qnj( z ) r _ Q Qn(r) j( z ) 

In addition, by ([75]) we know that for each fixed < r < m — 1 , 



lim 



Qn(r+l),ji Z ) \ ^ir^^p\ r ^ 

= F>'(z), z € K. CC\supp(aj) , 



nGA Q n (r),j(^) 

where I — (Ji; I2) is precisely the multi-index satisfying l\ = r mod (mi + 1), I2 = r mod (m.2 + 1). 
Therefore ([84]) follows. Relations ([85]) and (|86[) are proved analogously in view of (|80[) and ([82]) . □ 



Now, we need to introduce some notations. For j € {-m2, . . . , mi — 1}, set 

1, if Aj is to the left of Aj + i , 
— 1, if Aj is to the right of Aj + i . 

For multi-indices I = (h; I2) such that l\ + I2 > 2 we define 



A 



1, if 



83-U 



if 



j > h + 2 , 



< 



If /1 + l 2 = 1 then 



A 



3,1 



Sj-iSj, if j e {-i 2 + i,...,h - 1}, 
■Jj-, if je{-h- i,-h}, 

1, if j<_/ 2 _2. 



1, if j>h + 2, 

5j_i, if je{h,h + i}, 

-Sj, if j G {-I2- 1,-/2}, 

1, if j<-h~2, 
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and for l\ — l 2 = 



*i,(0;0) 



:= < 



1, 


if 


i>2, 


So, 


if 


.7 = 1: 


1, 


if 


J = 0, 


-S-u 


if 


J = -1 


1: 


if 


J < -2 



Recall that e nj denotes the sign of the varying measure 



dp n ,j(x) 



h nd (x)d<jj(x) 



Qn,j-l(x)Q a ,j+l(x) ' 

Lemma 11. For any n,n l G Z™ 1+1 (») x Z™ 2+1 (») and -m 2 < j < m 1 



(87) 



mi 

n Afe .i 

fc=J 



Proof. We will denote by sign(/, A) the sign of a function / on the interval A. Thus 

(88) 



&n l ,j . i Hn l .j Qnj'-l Qn.j + 1 . 

— - = sign( — — ^ — — , Ai 



-n,3 



^n,j Qn',j— 1 Qn ! ,j'+1 



If — Z 2 < j — 1 < Zi, then deg(<5 n i j^i) = 1 + deg(Q nj -_i) and, therefore, 
(89) 



sig n(i^zl, A.) = 



In 1 , 3— I 

If j — 1 < — Z 2 or j — 1 > ii, then deg(Q n ; = deg(Q nj _i), implying that 



(90) 



sign(- 



A,) = 1 



Analogously, we have that for — Z 2 < j + 1 < h 
(91) 

and for j + 1 < — / 2 or j + 1 > l± 



sign(^±l,A,) = -^ 



(92) 

From l(89] l -(|9T1 l it follows that 
(93) 



sign( 



Qn l 



,A,) = 1. 



/ Qn.j — l Qa,j+X a \ A 

Slgn( ?) — n — ' ^ = ' ' ' 



Now, by ([51]) 

Therefore, 
(94) 

Since 7i n ; mi 



r g^Lj+iW "n',j+l(') to i+lW 

~H a i,j(x) _ J x-t Q„i, J (*)Q n i, J+2 (*) 
J x-t Q»,j(t)Q n ,3+j(t) 



sigD.(n n ij/Hn,j,A 3j 



L = H n , mi = 1, the right hand side of fM)) equals 1 for j = m\. Hence ([57]) follows 
from (J55J), and (JSJ). □ 
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This lemma shows that E n ij/e n j depends on j, I, and the relative positions of the intervals Aj 
but not on n. Define the functions 

k=j+l \ K k I 

(the product should be understood to be equal to 1 when j = mi). 

Theorem 7. Let S 1 = N'{al, . . . , ff mi ), S 2 = Af'(o%, . . . , cr^J be given, and let A C Z" ll+1 (») x 
Z™ 2+1 («) be an infinite sequence of distinct multi-indices such that 

sup((m 2 + 1 )n 2 ,o — |n 2 |) < oo , sup((mi + 1)711,0 — |rii|) < oo . 

neA n£A 

Let us assume that there exists I = (Zij^jO < li < mi,0 < I2 < 7712, suc/i that for all n G A we 
Ziaue £/ia£ n' € x Z™ 2+1 (»). Let {-4 nj }™!_ m2 _i, n ^ ^> ^ e ^ e associated sequences of 

"monic" linear forms. Then, for each fixed j = — m 2 — f , . . . , mi, 

(95) lim = Af , /CcC\(supp(a,)Usup P ( ( 7 J+ i)) 

neA ^Anj^) J 
(supp(cr_ m2 _i) = SUpp(cr mi + i) = 0). 

Proof. It follows from definition of TL-aj and T~i n ',j that 

A n ij{z) _ £ n ' ,j+l n n> ,j ( z ) £nj + l + 1 Qn l ,j( z ) Qnj+l{z) 

•Aa,j( z ) e n j + ih n j(z) £ n ij + i K a l t j +1 Qn,j{z) Qn l ,j + l{ z ) 

By Lemma [FJ 

^±A^ = 1 1 K c c N supp( ) . 

Using Lemma [11] and Corollary [21 we have 

lim gn ' j+1 n ' i+1 = TT 

Finally, applying (175)) and (179p one obtains 

Q,»^)g^- + i(z) = w ^ c c N ( ( } u ( }) . 

neA Q n ,j{z) Qn',j+l{Z) J 

Putting these relations together we get (|55l) . □ 

Corollary 4. Let S 1 = Af'(<J^ ■ ■ ■ , <4J, S 2 = -^'(og, . . . ,o£j 6e groen, and /ei A C Z™ 1+1 (.) x 

Z? 2+ (•) 6e an infinite sequence of distinct multi-indices such that 

(96) sup((ra 2 + l)n2,o - |n 2 |) < 00 , sup((mi + l)ni, - |rii|) < 00 . 

neA neA 

Then, for each fixed j € {— m 2 , . . . , mi}, we have 

(97) lim " 4 " +p '; (z) = Af{z) , KCC\ (suppfo) U supp(a i+ i)) 

neA A-n : j(Z) ^ L 
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(supp(<7_ m2 _i) = supp((J TOl +i) = 0). Consequently, 



(98) 



lim |^l nj (z)| 1/|ni1 = TT ^(z)^ , K C CUsuppfaOUsuppfa+i)), 



2GL 



where m — lcm(mi + 1, m-i + 1). 

Proof. Using the same arguments employed to prove Corollary [3j we obtain (|97|) . From (|9T[) it is 
easy to deduce the |ni|-th root asymptotic of the linear forms. 

In fact, it is easy to see that for each n 6 A there exists no 6 Z™ 1+1 (») x Z™ 2+1 (») (which may 
depend on n), whose entries are uniformly bounded by a constant C independent of n (condition 
(|96|) is used), such that n = rp + no for some r S Z + . Write 

J / \ _ -^n,j( z ) A n - p j(z) -^Mo+pj ( z ) . / s 



•^n— p,j( z ) -^n— 2p,j ( z ) -^.rin , j (^) 



Then 



Obviously, 



1 1 1 r_1 

— log |^ n ,j(z)| - ,— -r log |-4ooj(«)| + |— T V log 



fe=0 



.4 



n + (fe+l 



v,j ( z ) 



.4 



no+fepj I 



1 



nSA | Hi 

and because of ([97 



lim r — f log \A a0l j(z)\ = 0, KcC\ (supp(cTj) U supp((T J+ i)), 



1 r_1 
lim 1 — r log 

nCA Til 



.4 



n + (fc+l) 



■An +kp,j {z) 



= — log 



, KcC\ (supp(cr i ) U supp(<r i+ i)), 



since |ni | = r|pi | + 0(1) — rm + 0(1), |ni | — > oo. 



□ 



The function appearing on the right hand side of (|9"8")l corresponds with the one on the right hand 
side of (|39[) associated to the vector equilibrium problem with interaction matrix C constructed 
taking pi^ = l/(mi + 1), < k < mi, and j»2,fe = 1/(^2 + 1), < k < m%. In that case, for each 
j = —rri2 — 1, ... , mi, we have 



G ^ 2 ) = II l<(*)l Vm . ^ e C \ (A,- U A J+1 ) 



(A- 



■m2-l — ^mi+l 



), where m = lcm(mi + 1, 7712 + 1). 



7. Application to mixed type Hermite Pade approximation 

Let S 1 = Af((jQ, ■ ■ ■ , a* ), S 2 = N{oq, . . . , of„ 2 ), 0q = a 2 be given. Let us introduce the row 
vectors 



U=(l,^ 1; ...,^ m2 ), V = 
and the (m2 + 1) x (mi + 1) dimensional matrix 

w = u'v, 



(1, s l x , . . . , s l mi ) 



42 U. FIDALGO, A. LOPEZ, G. LOPEZ, AND V. N. SOROKIN 

where the super-index t means taking transpose. Define the matrix Markov type function 

W(x)da%(x) 



§(z) _ 

z — X 

understanding that integration is carried out entry by entry on the matrix W. 

Fix ni = (m,o, ni,i,..., ni, mi ) € Z™ 1+1 and n 2 = (n 2 ,o, n 2 ,i, • . . , n 2 , m2 ) S Z™ 2+1 ,|n 2 | = 
| ni| — 1- It is easy to see that there exists a non zero vector polynomial 

A n = (a n ,o, ...,a„, mi ), deg(a„ )fe ) < ni, fe - 1, fc = 0, ...,mi, 

such that 

(99) §(z)A n (*) - B n (z) - (0(l/z" 2 ' o+1 ), . . . ,0(l/z n2 ^ +1 )) t =: <D(l/z n2+1 ), z -> oo, 

where D n — (rf n ,0i ■ ■ ■ , ^n,m 2 ) is some vector polynomial. When m 2 = 0, this construction is called 
type I Hermite-Pade approximation. If m\ = it is called of type II. When m± = m 2 = it 
reduces to diagonal Pade approximation. This definition is of mixed type. 

Lemma 12. For j = 0, . . . , m 2 , (s\ = 1) 

„ mi 

(100) j x v Y J ^,k{x)sl k {x)ds 2 j {x) = Q, v = 0,...,n 2tj -l. 

Proof. In fact, notice that according to (|99|). for each v 1 < v < n 2 j — 1, j = 0, . . . , m 2 , 

z " ( Jl anAz) J z-x ~ dnj(z) J = ° ( 1/Z ) ' Z ^ °°' 

(sf o = l) ano - the function on the left hand side is holomorphic in C\ Co(supp(<7o)). Using Lemma 
El we obtain (HOQl). □ 

Because of this Lemma, we see that A n is an n-th mixed type multiple orthogonal polynomial 
with respect to the pair (S 1 , S 2 ) and in the sequel we assume that it is "monic" . If 

in = (i»n,0 1 -"A,m s )> deg(6 nJ ) < n 2 ,j r- 1, j = 0, . . . , TO 2 , 

denotes a generic vector polynomial with the indicated degrees, (I100[) may be rewritten in matrix 
form as 

(101) J M n (x)W(x)A t n (x)dcrl(x) = 0, for all B n . 
Fix j G {0, . . . , m 2 }. For each fc 6 { — 1, . . . , — j — 1} define 

"i-{^C\ U^Ai : t/£ (z) < Uf(z), i = —1, . , . ,—j — l,i ^ fc}, n° a = C \ (A U A_i). 
Set 

Xj(z) := min{C/^(z) : fc = -1, . . . , -j - 1} 

and 

(Kn,o, K n , m2 y := S(z)A n (z) - O n (z). 
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Theorem 8. Let A = A(pi j0 , . . . ,Pi, mi ;P2 fi , ■ ■ ■ ,P2 >ma ) C Z™ 1+1 (.) x Z™ 2+1 (.), (S 1 ,^ 2 ) G Reg, 
■S 1 = AF{(Tq, . . . , ct^j ), and 5 2 = A/"(cr 2 , . . . , er 2 ^) &e given. Then for each j G {0, . . . , 7712} 



Urn l^-^l 1 /!^! = exp(- Xi (^)), £ C U^lM, 



(102) 
and 

(103) lira \HnA»)\ 1/]ni] < eM-Xi(z)), K. C C \ (U^ 1 A fe ). 
In particular, if P2,o = • • • = p^,m 2 = V( m 2 + 1); then 

(104) lim |H»,,(z)| 1/|ni1 = exp(-^ 1 (z)), KcC\ (U^ l A fc ). 

/I = /J(C) = (]I_ m , . . . ,~p m ) is t/ie equilibrium vector measure and (w^ m2 , . . . , a;^ ) is the system 
of equilibrium constants for the vector potential problem determined by the interaction matrix C 
defined in @ on the system of compact sets E}. = supp(cr^), k = 0, . . . , mi, Ek = supp(cr_ fe ), = 
-m 2 , ... ,0. 



Proof. Notice that (|101[) implies that 
g ( v A t (z) /" W(x)(A^(z) - A^z))^*) y WWA^j^fx) 

and taking 



/ 



C(l/z n2+i ), z^oo, 



>n(*) 



(^(A^-A^x))^) 



we obtain an integral expression for the remainder in 
Then 

(ftn,o(z), ■ ■ • ,^n,m 2 (^))* = 

In scalar form this says that 



W(x)A^(x)dcr 2 (a;) 



z — x 



-4n,o(x) , 2 



<2s,-(a;), i = 0, . ..,m 2 . 



Notice that (see (1721 ) 

^n.o(^) = A a -i(z). 

Let us establish a connection between the remainders lZ n j(z) and the forms ^4. n ,fc(^) with negative 
indices k G {—1, • ■ ■ , — j — 1}. 



Fix j G {1, . . . , m 2 }. We have 
and 

An.-j-l(z) - 

Consequently, 

(-iyn aij (z) - A^-j-xiz) 



An, ( x o) da%(x ) ■ ■ ■ da?(xj) 
(z - xq){x\ - x ) ■ ■ ■ (xj - Xj-i) ' 

A n ,o( x o) dcr%(x ) ■ ■ ■ d(T?(Xj) 
(x% - Xo) ■ ■ ■ {Xj - Xj-l)(z ~ Xj) 



-(xj - x Q )A n ,o(xo) da^(x Q ) ■ ■ -da 2 {xj) 

(z - X )(Xi - Xq) ■ ■ ■ (Xj - Xj-i)(z - Xj) ' 
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Since Xj — xq — xj — Xj-i + xj-i — ■ ■ ■ — x\ + x\ — xq, substituting this in the previous formula, 
we obtain 

k=0 

We have a triangular scheme of linear equations whose coefficients do not depend on n. We can 
solve for lZ n j in terms of A a -i, . . . , A n ^j-\. Using P5| one obtains that for each j G {0, . . . , 1712} 
(when j = the sum below is empty) 

n atj (z) = ^>i) fc -VL . _ a ^ z)An ^ k {z) + (-iyA a ,-j-i(z). 
fc=i 

Taking (|39[) into consideration, on Sl 3 _ k the term containing An.-k dominates the sum (notice that 
(a 2 ,, . . . , a 2 } (z) ^ 0, z G C \ A_fc) and (|102| immediately follows . On the complement of U^ili^ 
there is no dominating term and all we can conclude from the previous equality is (|103|) . 

Let p2,a = ■ ■ ■ = P2.m 2 = l/( m 2 + 1)- In this case, on C \ U^^Afc we have that /7^ 1 (z) < 
U%(z) <■■■ < U~ p j _ 1 (z) (see third sentence before Corollary[l|) and (fTM)) follows from (fTT]^|) . □ 

Remark 2. Fix j G {0, . . . , 7712}. For each k G { — 1, • • ■ , — j — 1} we could have defined 

Sl{ = {z G C \ UTf- 1 ^ : C/f (z) < Uf(z),i = -1, . . . , -j - 1, i± k}, n°_ x = C \ (So U £U). 

Taking into account that the polynomials Q ni i and the forms -4 n ,i may have at most one zero in 
each of the connected components of Aj \ Ei, one can prove in place of (|102|) - (|104[) convergence in 
capacity on each compact subset of the corresponding regions. □ 

We say that I\ C Z™ 1+1 (») is a complete, ordered, sequence of multi-indices if: 

a) For each n G Z + , there exists a unique ni G X\ such that |ni| = n. 

b) Any two multi-indices in X\ are ordered in the sense that all components of one of them are 
less than or equal to the corresponding components of the other one, or they are identical. 

Fix Xi C Z" ll+1 (»),X 2 C Z™ 2+1 (»), two complete, ordered sequences of multi-indices. Each 
n G Z + determines a unique ni G X\ and 112 G I2 such that n = |ni| ~ 1 112 1 + 1 • The corresponding 
"monic" mixed type multiple orthogonal polynomials we denote by A n . We can interchange the 
roles of the Nikishin systems S 1 ^ 2 , and determine a sequence of "monic" mixed type multiple 
orthogonal polynomials which we denote B„. It is easy to verify that the sequences {A„}, {B„}, n G 
Z + are bi-orthogonal. That is, 



(105) J n n ,(x)W(x)A t n (x)da 2 (x) 




The inequality in Q105P is a consequence of Lemma [5] With the same hypothesis, all the results of 
this paper hold true for the sequence {B„}, n G Z + . 
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